0.1 Load all of the throw data

# load omnibus dataframe
omnibus_df <- read_delim("../data/processed/omnibus/omnibus_raw.csv",
  delim = ",",
  col_types = cols(
    .default = col_double(),
    type = col_factor(),
    ppid = col_factor(),
    exp_label = col_factor(),
    experiment = col_factor(),
    hand = col_factor(),
    camera_tilt = col_factor(),
    surface_tilt = col_factor(),
    target = col_factor(),
    test_type = col_factor(),
    prior_anim = col_factor(),
    baseline_block = col_factor(),
    task_type = col_factor(),
    surface = col_factor(),
    anim_type = col_factor()
  )
) %>% # filter out practice blocks
  filter(block_num > 4)

# Optionally make learning rate summaries
# do the following if learning_rate_df.csv doesn't exist in ../data/processed
# This takes a loong time
if (!file.exists("../data/processed/learning_rate_df.csv")) {
  print(Sys.time())
  apply_exponentialFit <- function(df) {
    df %>%
      summarise(
        ppid = first(ppid),
        experiment = first(experiment),
        test_type = first(test_type),
        exponentialFit = exponentialFit2(norm_throw_deviation, mode = test_type[1])
      )
    }
  
  block_init_learning_rates <- omnibus_df %>%
    filter(str_detect(test_type, "init")) %>%
    group_by(ppid, experiment, test_type) %>%
    group_split() %>%
    future_map(apply_exponentialFit) %>%
    bind_rows() %>%
    unnest(cols = c('exponentialFit'))
  
  print("done")
  print(Sys.time())

  write_csv(block_init_learning_rates, "../data/processed/learning_rate_df.csv")
} else {
  print("learning_rate_df.csv exists, loading from file")
  block_init_learning_rates <- read_csv("../data/processed/learning_rate_df.csv",
                                        col_types = cols(
                                          .default = col_double(), 
                                          experiment = col_factor(),
                                          test_type = col_factor()
                                          ))
}
[1] "2023-06-23 19:18:08 EDT"
[1] "done"
[1] "2023-06-23 19:20:27 EDT"

1 Visualizing data (univariate)

Vectors representing the throw velocity (trace 0) and the velocity applied to the ball (trace 1). The y dimention of the throw is essentially ignored (in reality there is a slight tilt added to account for the tilt of the surface).

test_ppt <- 3

test_df <- omnibus_df %>% filter(ppid == test_ppt)


trial <- 250
trial_df <- filter(test_df, trial_num == trial)

x <- trial_df$flick_velocity_x
y <- trial_df$flick_velocity_y
z <- trial_df$flick_velocity_z

x2 <- trial_df$flick_direction_x * -1
y2 <- trial_df$flick_direction_y * -1
z2 <- trial_df$flick_direction_z * -1

# plot both
p <- plot_ly(x = c(0, x), y = c(0, y), z = c(0, z), type = "scatter3d", mode = "lines") %>%
  add_trace(x = c(0, x2), y = c(0, y2), z = c(0, z2), type = "scatter3d", mode = "lines") %>%
  layout(scene = list(
    xaxis = list(title = "x", range = c(-2, 2)),
    yaxis = list(title = "y", range = c(-1, 3)),
    zaxis = list(title = "z", range = c(-1, 3))
  ))

# Render the plot
p

note: this is a rotated trial

1.0.1 Distribution of errors

# plot distribution of error_size
p <- ggplot(omnibus_df, aes(
  x = error_size,
  fill = type
)) +
  geom_histogram(binwidth = .5, alpha = .6) +
  theme_minimal() +
  theme(text = element_text(size = 11)) +
  scale_fill_manual(values = c("#f9982c", "#d40000")) +
  labs(x = "Error Size (cm)", y = "Count")

p

1.0.2 Distribution of throw angles

# plot distribution of error_size
p <- ggplot(omnibus_df, aes(
  x = throw_deviation,
  fill = type
)) +
  geom_histogram(binwidth = 1, alpha = .6) +
  theme_minimal() +
  theme(text = element_text(size = 11)) +
  scale_fill_manual(values = c("#f9982c", "#d40000")) +
  labs(
    x = "Throw Angle (°)", y = "Count"
  ) + # dashed lines at 0, -15, -30
  geom_vline(
    xintercept = c(0, -15, -30), linewidth = 0.4,
    colour = "#CCCCCC", linetype = "dashed"
  ) + # ticks of 15 degrees
  scale_x_continuous(
    breaks = c(-30, 0, 30, -60, -90)
  )

p

2 Original Experiments

2.0.1 Plot ANGULAR DEVIATIONS (hand angles)

Note: Blues = Acceleration Perturbations

# rest of the exps
data_per_group <- omnibus_df %>%
  filter(exp_label == "original_exps" | exp_label == "curved_path") %>%
  group_by(experiment, test_type, trial_num) %>%
  summarise(
    mean_deviation = mean(throw_deviation),
    ci_deviation = vector_confint(throw_deviation),
    .groups = "drop"
  )

# set up plot
p <- data_per_group %>%
  ggplot(
    aes(
      x = trial_num, y = mean_deviation, colour = experiment
    )
  ) +
  theme_classic() +
  # theme(legend.position = "none") +
  labs(
    x = "Trial Number",
    y = "Throw Angle (°)"
  )

# add horizontal lines
p <- p +
  geom_hline(
    yintercept = c(0, -30), linewidth = 0.4,
    colour = "#CCCCCC", linetype = "solid"
  ) +
  geom_hline(
    yintercept = c(-15), linewidth = 0.4,
    colour = "#CCCCCC", linetype = "dashed"
  )

# p <- p +
#   scale_y_continuous(
#     limits = c(-10, 35),
#     breaks = c(0, 15, 30),
#     labels = c(0, 15, 30)
#   ) +
#   scale_x_continuous(
#     limits = c(0, 180),
#     breaks = c(0, 60, 120, 180),
#     labels = c(0, 60, 120, 180)
#   )

# set font size to 11
p <- p +
  theme(text = element_text(size = 11))

# add confidence intervals and data points
p <- p + geom_ribbon(
  aes(
    ymin = mean_deviation - ci_deviation,
    ymax = mean_deviation + ci_deviation,
    fill = experiment
  ),
  colour = NA, alpha = 0.3
) + geom_line()

# set colour palette
p <- p +
  scale_colour_manual(values = pallete_list) +
  scale_fill_manual(values = pallete_list)


# # save
# if (save_plots) {
#   ggsave(
#   p,
#   filename = "../plots/paper_figs/sr_30_training.pdf", device = "pdf",
#   height = 4, width = 6
#   )
#   }

ggplotly(p)

# p

2.0.2 Trial sets of interest only

Note: Blues = Acceleration Perturbations

# filter out just the trials of interest
data_per_group <- data_per_group %>%
  filter(
    test_type != "other"
  )
# add a dummy column with repeating sequence
# NOTE: this can't be combined with above since we are using nrow
data_per_group <- data_per_group %>%
  mutate(dummy_x = rep(1:(nrow(data_per_group) / NUM_EXPS),
    length.out = nrow(data_per_group)
  ))

# set up plot
p <- data_per_group %>%
  ggplot(
    aes(
      x = dummy_x, y = mean_deviation, colour = experiment
    )
  ) +
  theme_classic() +
  # theme(legend.position = "none") +
  labs(
    x = "Trial Number",
    y = "Throw Angle (°)"
  )

# add horizontal lines
p <- p +
  geom_hline(
    yintercept = c(0, -30), linewidth = 0.4,
    colour = "#CCCCCC", linetype = "solid"
  ) +
  geom_hline(
    yintercept = c(-15), linewidth = 0.4,
    colour = "#CCCCCC", linetype = "dashed"
  )

# add confidence intervals and data points
for (unique_test_type in unique(data_per_group$test_type)) {
  # get the data for this block
  to_plot_data <- filter(data_per_group, test_type == unique_test_type)

  p <- p + geom_ribbon(
    data = to_plot_data,
    aes(
      ymin = mean_deviation - ci_deviation,
      ymax = mean_deviation + ci_deviation,
      fill = experiment
    ), colour = NA, alpha = 0.3
  ) + geom_line(
    data = to_plot_data
  )
}

# set colour palette
p <- p +
  scale_colour_manual(values = pallete_list) +
  scale_fill_manual(values = pallete_list)

ggplotly(p)

2.0.3 Plot NORMALIZED ANGULAR DEVIATIONS (hand angles)

Note: Blues = Acceleration Perturbations

# rest of the exps
data_per_group <- omnibus_df %>%
  filter(exp_label == "original_exps" | exp_label == "curved_path") %>%
  group_by(experiment, test_type, trial_num) %>%
  summarise(
    mean_deviation = mean(norm_throw_deviation),
    ci_deviation = vector_confint(norm_throw_deviation),
    .groups = "drop"
  )

# set up plot
p <- data_per_group %>%
  ggplot(
    aes(
      x = trial_num, y = mean_deviation, colour = experiment
    )
  ) +
  theme_classic() +
  # theme(legend.position = "none") +
  labs(
    x = "Trial Number",
    y = "Normalized Throw Angle"
  )

# add horizontal lines
p <- p +
  geom_hline(
    yintercept = c(0, 1, 2), linewidth = 0.4,
    colour = "#CCCCCC", linetype = "solid"
  ) +
  geom_hline(
    yintercept = c(0.5, 1.5), linewidth = 0.4,
    colour = "#CCCCCC", linetype = "dashed"
  )

# p <- p +
#   scale_y_continuous(
#     limits = c(-10, 35),
#     breaks = c(0, 15, 30),
#     labels = c(0, 15, 30)
#   ) +
#   scale_x_continuous(
#     limits = c(0, 180),
#     breaks = c(0, 60, 120, 180),
#     labels = c(0, 60, 120, 180)
#   )

# set font size to 11
p <- p +
  theme(text = element_text(size = 11))

# add confidence intervals and data points
p <- p + geom_ribbon(
  aes(
    ymin = mean_deviation - ci_deviation,
    ymax = mean_deviation + ci_deviation,
    fill = experiment
  ),
  colour = NA, alpha = 0.3
) + geom_line()

# set colour palette
p <- p +
  scale_colour_manual(values = pallete_list) +
  scale_fill_manual(values = pallete_list)


# # save
# if (save_plots) {
#   ggsave(
#   p,
#   filename = "../plots/paper_figs/sr_30_training.pdf", device = "pdf",
#   height = 4, width = 6
#   )
#   }

ggplotly(p)

# p

2.0.4 Trial sets of interest only

Note: Blues = Acceleration Perturbations

# filter out just the trials of interest
data_per_group <- data_per_group %>%
  filter(
    test_type != "other"
  )
# add a dummy column with repeating sequence
# NOTE: this can't be combined with above since we are using nrow
data_per_group <- data_per_group %>%
  mutate(dummy_x = rep(1:(nrow(data_per_group) / NUM_EXPS),
    length.out = nrow(data_per_group)
  ))

# set up plot
p <- data_per_group %>%
  ggplot(
    aes(
      x = dummy_x, y = mean_deviation, colour = experiment
    )
  ) +
  theme_classic() +
  # theme(legend.position = "none") +
  labs(
    x = "Trial Number",
    y = "Normalized Throw Angle"
  )

# add horizontal lines
p <- p +
  geom_hline(
    yintercept = c(0, 1, 2), linewidth = 0.4,
    colour = "#CCCCCC", linetype = "solid"
  ) +
  geom_hline(
    yintercept = c(0.5, 1.5), linewidth = 0.4,
    colour = "#CCCCCC", linetype = "dashed"
  )

# add confidence intervals and data points
for (unique_test_type in unique(data_per_group$test_type)) {
  # get the data for this block
  to_plot_data <- filter(data_per_group, test_type == unique_test_type)

  p <- p + geom_ribbon(
    data = to_plot_data,
    aes(
      ymin = mean_deviation - ci_deviation,
      ymax = mean_deviation + ci_deviation,
      fill = experiment
    ), colour = NA, alpha = 0.3
  ) + geom_line(
    data = to_plot_data
  )
}

# set colour palette
p <- p +
  scale_colour_manual(values = pallete_list) +
  scale_fill_manual(values = pallete_list)

ggplotly(p)

2.1 Learning Rates, Start Points and Asymptotes

2.1.1 Learning Rates

data_ppt <- block_init_learning_rates %>%
  filter(experiment != "a_ball_roll_animate_surface")

data_group <- data_ppt %>%
  group_by(experiment, test_type) %>%
  summarise(
    mean_learning_rate = mean(exp_fit_lambda),
    ci_learning_rate = vector_confint(exp_fit_lambda),
    mean_high = mean(exp_fit_N0),
    ci_high = vector_confint(exp_fit_N0),
    mean_low = mean(exp_fit_displace),
    ci_low = vector_confint(exp_fit_displace),
    .groups = "drop"
  )

p <- data_group %>%
  ggplot(
    aes(x = experiment, y = mean_learning_rate, colour = experiment)
  ) +
  theme_classic() +
  labs(
    x = NULL,
    y = "Learning Rate"
  ) +
  facet_wrap(~test_type)

# remove all x axis labels
p <- p + theme(axis.text.x = element_blank())

# for the colour legend, only show the first 7 Note this doesn't work for the plotly plot
p <- p + guides(colour = guide_legend(override.aes = list(alpha = 1)))

# add data points
p <- p +
  geom_beeswarm(
    data = data_ppt,
    aes(
      y = exp_fit_lambda
    ),
    alpha = 0.1,
    size = 1
  ) +
  geom_linerange(aes(
    ymin = mean_learning_rate - ci_learning_rate,
    ymax = mean_learning_rate + ci_learning_rate
  ), alpha = 0.5, lwd = 2) +
  geom_point()

# set colour palette
p <- p +
  scale_colour_manual(values = pallete_list) +
  scale_fill_manual(values = pallete_list)

ggplotly(p)
# p

2.1.2 High Points (Asymptotes or Starts)

p <- data_group %>%
  ggplot(
    aes(x = experiment, y = mean_high, colour = experiment)
  ) +
  theme_classic() +
  labs(
    x = NULL,
    y = "High Point"
  ) +
  facet_wrap(~test_type)

# remove all x axis labels
p <- p + theme(axis.text.x = element_blank())

# for the colour legend, only show the first 7 Note this doesn't work for the plotly plot
p <- p + guides(colour = guide_legend(override.aes = list(alpha = 1)))

# add data points
p <- p +
  geom_beeswarm(
    data = data_ppt,
    aes(
      y = exp_fit_N0
    ),
    alpha = 0.1,
    size = 1
  ) +
  geom_linerange(aes(
    ymin = mean_high - ci_high,
    ymax = mean_high + ci_high
  ), alpha = 0.5, lwd = 2) +
  geom_point()

# set colour palette
p <- p +
  scale_colour_manual(values = pallete_list) +
  scale_fill_manual(values = pallete_list)

ggplotly(p)
# p

2.1.3 Low Points (Asymptotes or Starts)


### Low Points (Asymptotes or Starts)

p <- data_group %>%
  ggplot(
    aes(x = experiment, y = mean_low, colour = experiment)
  ) +
  theme_classic() +
  labs(
    x = NULL,
    y = "Low Point"
  ) +
  facet_wrap(~test_type)

# remove all x axis labels
p <- p + theme(axis.text.x = element_blank())

# for the colour legend, only show the first 7 Note this doesn't work for the plotly plot
p <- p + guides(colour = guide_legend(override.aes = list(alpha = 1)))

# add data points
p <- p +
  geom_beeswarm(
    data = data_ppt,
    aes(
      y = exp_fit_displace
    ),
    alpha = 0.1,
    size = 1
  ) +
  geom_linerange(aes(
    ymin = mean_low - ci_low,
    ymax = mean_low + ci_low
  ), alpha = 0.5, lwd = 2) +
  geom_point()

# set colour palette
p <- p +
  scale_colour_manual(values = pallete_list) +
  scale_fill_manual(values = pallete_list)

ggplotly(p)
# p

For washout: The CUED accel + curved have a lower starting point (therefore – cue works). VMR group has slightly lower. When comparing everything with a high starting point, the ACCEL group has a much FASTER learning rate. When transferring, no difference in learning rates.

# Compare learning rates and high points
p <- data_ppt %>% 
  ggplot(
    aes(x = exp_fit_lambda, y = exp_fit_N0, colour = experiment)
  ) +
  theme_classic() +
  labs(
    x = "Learning Rate",
    y = "High Point"
  ) +
  facet_wrap(~test_type) +
  geom_point() +
  scale_colour_manual(values = pallete_list) +
  scale_fill_manual(values = pallete_list)

ggplotly(p)

2.1.4 Correlations within the 3 variables above


# 3d plot of learning rates, high points, and low points
# filter out just the washout trials for now
data_ppt <- block_init_learning_rates %>%
  filter(test_type == "washout_init")

p <- plot_ly(type = 'scatter3d', mode = 'markers')

for (experiment_ in unique(data_ppt$experiment)) {
  # get the data for this block
  to_plot_data <- filter(data_ppt, experiment == experiment_)
  
  p <- p %>%
    add_trace(
      x = to_plot_data$exp_fit_lambda, 
      y = to_plot_data$exp_fit_N0, 
      z = to_plot_data$exp_fit_displace, 
      color = I(pallete_list[experiment_]),
      name = experiment_
    )
}

p
  • We might want to test if data points cluster differently depending on the experiment conditions.
  • Note: It is a good idea to do some EDA to check if all 3 dimensions are required. Correlation matrix?
### TESTING
# rest of the exps
data_ <- omnibus_df %>%
  filter(experiment == "accel_uncued", test_type == "washout_init")
  
# set up plot
p <- data_ %>%
  ggplot(
    aes(
      x = trial_num, y = norm_throw_deviation, colour = ppid
    )
  ) +
  theme_classic() +
  # theme(legend.position = "none") +
  labs(
    x = "Trial Number",
    y = "Normalized Throw Angle"
  )

# add horizontal lines
p <- p +
  geom_hline(
    yintercept = c(0, 1, 2), linewidth = 0.4,
    colour = "#CCCCCC", linetype = "solid"
  ) +
  geom_hline(
    yintercept = c(0.5, 1.5), linewidth = 0.4,
    colour = "#CCCCCC", linetype = "dashed"
  )

# p <- p +
#   scale_y_continuous(
#     limits = c(-10, 35),
#     breaks = c(0, 15, 30),
#     labels = c(0, 15, 30)
#   ) +
#   scale_x_continuous(
#     limits = c(0, 180),
#     breaks = c(0, 60, 120, 180),
#     labels = c(0, 60, 120, 180)
#   )

# set font size to 11
p <- p +
  theme(text = element_text(size = 11))

# add confidence intervals and data points
p <- p + geom_line()



# # save
# if (save_plots) {
#   ggsave(
#   p,
#   filename = "../plots/paper_figs/sr_30_training.pdf", device = "pdf",
#   height = 4, width = 6
#   )
#   }

ggplotly(p)

# p

2.2 Plot ERROR SIZE

Note: Blues = Acceleration Perturbations

# original experiments only
data_per_group <- omnibus_df %>%
  filter(exp_label == "original_exps" | exp_label == "curved_path") %>%
  group_by(experiment, test_type, trial_num) %>%
  summarise(
    mean_deviation = mean(error_size),
    ci_deviation = vector_confint(error_size),
    .groups = "drop"
  )

# set up plot
p <- data_per_group %>%
  ggplot(
    aes(
      x = trial_num, y = mean_deviation, colour = experiment
    )
  ) +
  theme_classic() +
  # theme(legend.position = "none") +
  labs(
    x = "Trial Number",
    y = "Absolute Target Error (cm)"
  )

# add horizontal lines
p <- p +
  geom_hline(
    yintercept = c(0, 40), linewidth = 0.4,
    colour = "#CCCCCC", linetype = "solid"
  ) +
  geom_hline(
    yintercept = c(20), linewidth = 0.4,
    colour = "#CCCCCC", linetype = "dashed"
  )

# p <- p +
#   scale_y_continuous(
#     limits = c(-10, 35),
#     breaks = c(0, 15, 30),
#     labels = c(0, 15, 30)
#   ) +
#   scale_x_continuous(
#     limits = c(0, 180),
#     breaks = c(0, 60, 120, 180),
#     labels = c(0, 60, 120, 180)
#   )

# set font size to 11
p <- p +
  theme(text = element_text(size = 11))

# add confidence intervals and data points
p <- p + geom_ribbon(
  aes(
    ymin = mean_deviation - ci_deviation,
    ymax = mean_deviation + ci_deviation,
    fill = experiment
  ),
  colour = NA, alpha = 0.3
) + geom_line()

# set colour palette
p <- p +
  scale_colour_manual(values = pallete_list) +
  scale_fill_manual(values = pallete_list)


# # save
# if (save_plots) {
#   ggsave(
#   p,
#   filename = "../plots/paper_figs/sr_30_training.pdf", device = "pdf",
#   height = 4, width = 6
#   )
#   }

ggplotly(p)

# p

visible vs non-visible tilt doesn’t affect the 15-degree rotation condition. But affects all other conditions. So 15-degree rotation

2.2.1 Trial sets of interest only

Note: Blues = Acceleration Perturbations

# filter out just the trials of interest
data_per_group <- data_per_group %>%
  filter(
    test_type != "other"
  )
# add a dummy column with repeating sequence
# NOTE: this can't be combined with above since we are using nrow
data_per_group <- data_per_group %>%
  mutate(dummy_x = rep(1:(nrow(data_per_group) / NUM_EXPS),
    length.out = nrow(data_per_group)
  ))

# set up plot
p <- data_per_group %>%
  ggplot(
    aes(
      x = dummy_x, y = mean_deviation, colour = experiment
    )
  ) +
  theme_classic() +
  # theme(legend.position = "none") +
  labs(
    x = "Trial Number",
    y = "Absolute Target Error (cm)"
  )

# add horizontal lines
p <- p +
  geom_hline(
    yintercept = c(0, 40), linewidth = 0.4,
    colour = "#CCCCCC", linetype = "solid"
  ) +
  geom_hline(
    yintercept = c(20), linewidth = 0.4,
    colour = "#CCCCCC", linetype = "dashed"
  )

# add confidence intervals and data points
for (unique_test_type in unique(data_per_group$test_type)) {
  # get the data for this block
  to_plot_data <- filter(data_per_group, test_type == unique_test_type)

  p <- p + geom_ribbon(
    data = to_plot_data,
    aes(
      ymin = mean_deviation - ci_deviation,
      ymax = mean_deviation + ci_deviation,
      fill = experiment
    ), colour = NA, alpha = 0.3
  ) + geom_line(
    data = to_plot_data
  )
}

# set colour palette
p <- p +
  scale_colour_manual(values = pallete_list) +
  scale_fill_manual(values = pallete_list)

ggplotly(p)

3 Animated Surface Follow-up

3.0.1 Plot ANGULAR DEVIATIONS (hand angles)

# isolate animate_surface exp
data_per_group <- omnibus_df %>%
  filter(exp_label == "animate_surface") %>%
  group_by(prior_anim, block_num, trial_num_in_block, trial_num) %>%
  summarise(
    mean_deviation = mean(throw_deviation),
    ci_deviation = vector_confint(throw_deviation)
  )

# order the factors for assigning colour pallets
data_per_group$prior_anim <- factor(
  data_per_group$prior_anim,
  levels = c(
    "none", "half_anim", "full_anim"
  )
)

# set up plot
p <- data_per_group %>%
  ggplot(
    aes(
      x = trial_num, y = mean_deviation,
      ymin = mean_deviation - ci_deviation,
      ymax = mean_deviation + ci_deviation
    )
  ) +
  theme_classic() +
  # theme(legend.position = "none") +
  labs(
    x = "Trial Number",
    y = "Throw Angle (°)"
  )

# add horizontal lines
p <- p +
  geom_hline(
    yintercept = c(0, -30), linewidth = 0.4,
    colour = "#CCCCCC", linetype = "solid"
  ) +
  geom_hline(
    yintercept = c(-15), linewidth = 0.4,
    colour = "#CCCCCC", linetype = "dashed"
  )

# p <- p +
#   scale_y_continuous(
#     limits = c(-10, 35),
#     breaks = c(0, 15, 30),
#     labels = c(0, 15, 30)
#   ) +
#   scale_x_continuous(
#     limits = c(0, 180),
#     breaks = c(0, 60, 120, 180),
#     labels = c(0, 60, 120, 180)
#   )

# set font size to 11
p <- p +
  theme(text = element_text(size = 11))

# repeat for prior_anim == "half", "full" and "wait"
for (unique_prior_anim in unique(data_per_group$prior_anim)) {
  # get the data for this block
  to_plot_data <- filter(data_per_group, prior_anim == unique_prior_anim)
  # loop through the unique blocks in to_plot_data
  for (block in unique(to_plot_data$block_num)) {
    # get the data for this block
    block_data <- filter(to_plot_data, block_num == block)
    # add the data, use the pallete_list to get the colour
    p <- p + geom_ribbon(
      data = block_data,
      aes(fill = prior_anim),
      colour = NA, alpha = 0.3
    ) + geom_line(
      data = block_data,
      aes(colour = prior_anim)
    )
  }
}

# set colour palette
p <- p +
  scale_colour_manual(values = pallete_list) +
  scale_fill_manual(values = pallete_list)

# # save
# if (save_plots) {
#   ggsave(
#   p,
#   filename = "../plots/paper_figs/sr_30_training.pdf", device = "pdf",
#   height = 4, width = 6
#   )
#   }

ggplotly(p)
# p

3.0.2 Washout trials after half VS full animations

# first, isolate the data
data <- omnibus_df %>%
  filter(
    exp_label == "animate_surface",
    baseline_block == FALSE,
    test_type == "washout_anim"
  )

data_per_ppt <- data %>%
  group_by(ppid, prior_anim, trial_num_in_block) %>%
  summarise(
    ppt_mean_deviation = median(throw_deviation),
    ppt_ci_deviation = vector_confint(throw_deviation),
    n = n()
  )

data_per_group <- data_per_ppt %>%
  group_by(prior_anim, trial_num_in_block) %>%
  summarise(
    mean_deviation = mean(ppt_mean_deviation),
    ci_deviation = vector_confint(ppt_mean_deviation),
    n = sum(n)
  )

# set up plot
p <- data_per_group %>%
  ggplot(
    aes(
      x = trial_num_in_block, y = mean_deviation,
      colour = prior_anim, fill = prior_anim
    )
  ) +
  theme_classic() +
  # theme(legend.position = "none") +
  labs(
    x = "Trial Number in Block",
    y = "Throw Angle (°)"
  )

# add horizontal lines
p <- p +
  geom_hline(
    yintercept = c(0, -30), linewidth = 0.4,
    colour = "#CCCCCC", linetype = "solid"
  ) +
  geom_hline(
    yintercept = c(-15), linewidth = 0.4,
    colour = "#CCCCCC", linetype = "dashed"
  )

# add data points
p <- p + geom_beeswarm(
  data = data_per_ppt,
  aes(
    y = ppt_mean_deviation,
    colour = prior_anim
  ),
  size = 1, dodge.width = 0.5
) + geom_ribbon(
  aes(
    ymin = mean_deviation - ci_deviation,
    ymax = mean_deviation + ci_deviation
  ),
  colour = NA, alpha = 0.3
) + geom_line()


# set colour palette
p <- p +
  scale_colour_manual(values = pallete_list) +
  scale_fill_manual(values = pallete_list)

ggplotly(p)

4 Deprecated (Don’t Run)

4.1 Success manifolds

4.1.1 Without any tilts

# ggplotly(plot_success_manifold_no_tilt())
plot_success_manifold_no_tilt()

4.1.2 With tilt present

ggplotly(plot_success_manifold_tilt())
LS0tDQp0aXRsZTogIkJpbGxpYXJkcyBhbmQgVGlsdHMgQW5hbHlzaXMgTm90ZWJvb2siDQphdXRob3I6ICJTaGFuYWF0aGFuYW4gTW9kY2hhbGluZ2FtIg0Kb3V0cHV0OiANCiAgaHRtbF9ub3RlYm9vazoNCiAgICB0b2M6IHRydWUNCiAgICB0b2NfZmxvYXQ6IHRydWUNCiAgICBudW1iZXJfc2VjdGlvbnM6IHRydWUNCiAgICBkZl9wcmludDogcGFnZWQNCiAgICBjb2RlX2ZvbGRpbmc6IGhpZGUNCi0tLQ0KDQpgYGB7ciBzZXR1cCwgaW5jbHVkZT1GQUxTRSwgd2FybmluZz1GQUxTRX0NCnJtKGxpc3QgPSBscygpKSAjIGNsZWFuIGVudmlyb25tZW50DQoNCnNvdXJjZSgiLi4vc3JjL2hlbHBlcl9mdW5jcy5SIikNCnNvdXJjZSgiLi4vc2NyaXB0cy9maWd1cmVfZnVuY3MuUiIpDQpsaWJyYXJ5KGRhdGEudGFibGUpDQpsaWJyYXJ5KHRpZHl2ZXJzZSkNCmxpYnJhcnkoZ2diZWVzd2FybSkNCmxpYnJhcnkoZXopICMgZm9yIEFOT1ZBcw0KbGlicmFyeShlZmZlY3RzaXplKSAjIGZvciBldGEtc3F1YXJlZA0KbGlicmFyeShwbG90bHkpDQoNCmxpYnJhcnkoZnVycnIpDQpwbGFuKG11bHRpc2Vzc2lvbikNCg0Kb3B0aW9ucyhkcGx5ci5zdW1tYXJpc2UuaW5mb3JtID0gRikNCg0KIyB2YXJzDQpvbW5pYnVzX3BhdGggPC0gIi4uL2RhdGEvcHJvY2Vzc2VkL29tbmlidXMvb21uaWJ1c19yYXcuY3N2Ig0KDQojIGNvbnZlcnQgdGhlIGFib3ZlIGludG8gYSBsaXN0DQpwYWxsZXRlX2xpc3QgPC0gYygNCiAgInJvdDMwX2N1ZWRfdGlsdCIgPSAiI2Q0MDAwMCIsDQogICJyb3QzMF91bmN1ZWQiID0gIiNmOTk4MmMiLA0KICAiYWNjZWxfY3VlZF90aWx0IiA9ICIjMDc1MDliIiwNCiAgImFjY2VsX3VuY3VlZCIgPSAiIzVmYjY5NiIsDQogICJjdXJ2ZWRfY3VlZF90aWx0IiA9ICIjMmI1NzQ3IiwNCiAgInJvdDE1X2N1ZWRfdGlsdCIgPSAiIzc3MDIwMiIsDQogICJyb3QxNV91bmN1ZWQiID0gIiNhNzYzMTUiLA0KICAibm9uZSIgPSAiI2Y5OTgyYyIsDQogICJoYWxmX2FuaW0iID0gIiM1ZmI2OTYiLA0KICAiZnVsbF9hbmltIiA9ICIjMDc1MDliIiwNCiAgIndhaXQiID0gIiNhNzYzMTUiDQopDQoNCk5VTV9FWFBTIDwtIDcNCmBgYA0KDQoNCg0KIyMgTG9hZCBhbGwgb2YgdGhlIHRocm93IGRhdGENCg0KDQoNCmBgYHtyfQ0KIyBsb2FkIG9tbmlidXMgZGF0YWZyYW1lDQpvbW5pYnVzX2RmIDwtIHJlYWRfZGVsaW0oIi4uL2RhdGEvcHJvY2Vzc2VkL29tbmlidXMvb21uaWJ1c19yYXcuY3N2IiwNCiAgZGVsaW0gPSAiLCIsDQogIGNvbF90eXBlcyA9IGNvbHMoDQogICAgLmRlZmF1bHQgPSBjb2xfZG91YmxlKCksDQogICAgdHlwZSA9IGNvbF9mYWN0b3IoKSwNCiAgICBwcGlkID0gY29sX2ZhY3RvcigpLA0KICAgIGV4cF9sYWJlbCA9IGNvbF9mYWN0b3IoKSwNCiAgICBleHBlcmltZW50ID0gY29sX2ZhY3RvcigpLA0KICAgIGhhbmQgPSBjb2xfZmFjdG9yKCksDQogICAgY2FtZXJhX3RpbHQgPSBjb2xfZmFjdG9yKCksDQogICAgc3VyZmFjZV90aWx0ID0gY29sX2ZhY3RvcigpLA0KICAgIHRhcmdldCA9IGNvbF9mYWN0b3IoKSwNCiAgICB0ZXN0X3R5cGUgPSBjb2xfZmFjdG9yKCksDQogICAgcHJpb3JfYW5pbSA9IGNvbF9mYWN0b3IoKSwNCiAgICBiYXNlbGluZV9ibG9jayA9IGNvbF9mYWN0b3IoKSwNCiAgICB0YXNrX3R5cGUgPSBjb2xfZmFjdG9yKCksDQogICAgc3VyZmFjZSA9IGNvbF9mYWN0b3IoKSwNCiAgICBhbmltX3R5cGUgPSBjb2xfZmFjdG9yKCkNCiAgKQ0KKSAlPiUgIyBmaWx0ZXIgb3V0IHByYWN0aWNlIGJsb2Nrcw0KICBmaWx0ZXIoYmxvY2tfbnVtID4gNCkNCg0KIyBPcHRpb25hbGx5IG1ha2UgbGVhcm5pbmcgcmF0ZSBzdW1tYXJpZXMNCiMgZG8gdGhlIGZvbGxvd2luZyBpZiBsZWFybmluZ19yYXRlX2RmLmNzdiBkb2Vzbid0IGV4aXN0IGluIC4uL2RhdGEvcHJvY2Vzc2VkDQojIFRoaXMgdGFrZXMgYSBsb29uZyB0aW1lDQppZiAoIWZpbGUuZXhpc3RzKCIuLi9kYXRhL3Byb2Nlc3NlZC9sZWFybmluZ19yYXRlX2RmLmNzdiIpKSB7DQogIHByaW50KFN5cy50aW1lKCkpDQogIGFwcGx5X2V4cG9uZW50aWFsRml0IDwtIGZ1bmN0aW9uKGRmKSB7DQogICAgZGYgJT4lDQogICAgICBzdW1tYXJpc2UoDQogICAgICAgIHBwaWQgPSBmaXJzdChwcGlkKSwNCiAgICAgICAgZXhwZXJpbWVudCA9IGZpcnN0KGV4cGVyaW1lbnQpLA0KICAgICAgICB0ZXN0X3R5cGUgPSBmaXJzdCh0ZXN0X3R5cGUpLA0KICAgICAgICBleHBvbmVudGlhbEZpdCA9IGV4cG9uZW50aWFsRml0Mihub3JtX3Rocm93X2RldmlhdGlvbiwgbW9kZSA9IHRlc3RfdHlwZVsxXSkNCiAgICAgICkNCiAgICB9DQogIA0KICBibG9ja19pbml0X2xlYXJuaW5nX3JhdGVzIDwtIG9tbmlidXNfZGYgJT4lDQogICAgZmlsdGVyKHN0cl9kZXRlY3QodGVzdF90eXBlLCAiaW5pdCIpKSAlPiUNCiAgICBncm91cF9ieShwcGlkLCBleHBlcmltZW50LCB0ZXN0X3R5cGUpICU+JQ0KICAgIGdyb3VwX3NwbGl0KCkgJT4lDQogICAgZnV0dXJlX21hcChhcHBseV9leHBvbmVudGlhbEZpdCkgJT4lDQogICAgYmluZF9yb3dzKCkgJT4lDQogICAgdW5uZXN0KGNvbHMgPSBjKCdleHBvbmVudGlhbEZpdCcpKQ0KICANCiAgcHJpbnQoImRvbmUiKQ0KICBwcmludChTeXMudGltZSgpKQ0KDQogIHdyaXRlX2NzdihibG9ja19pbml0X2xlYXJuaW5nX3JhdGVzLCAiLi4vZGF0YS9wcm9jZXNzZWQvbGVhcm5pbmdfcmF0ZV9kZi5jc3YiKQ0KfSBlbHNlIHsNCiAgcHJpbnQoImxlYXJuaW5nX3JhdGVfZGYuY3N2IGV4aXN0cywgbG9hZGluZyBmcm9tIGZpbGUiKQ0KICBibG9ja19pbml0X2xlYXJuaW5nX3JhdGVzIDwtIHJlYWRfY3N2KCIuLi9kYXRhL3Byb2Nlc3NlZC9sZWFybmluZ19yYXRlX2RmLmNzdiIsDQogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgY29sX3R5cGVzID0gY29scygNCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIC5kZWZhdWx0ID0gY29sX2RvdWJsZSgpLCANCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIGV4cGVyaW1lbnQgPSBjb2xfZmFjdG9yKCksDQogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICB0ZXN0X3R5cGUgPSBjb2xfZmFjdG9yKCkNCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICkpDQp9DQpgYGANCg0KIyBWaXN1YWxpemluZyBkYXRhICh1bml2YXJpYXRlKQ0KVmVjdG9ycyByZXByZXNlbnRpbmcgdGhlIHRocm93IHZlbG9jaXR5ICh0cmFjZSAwKSBhbmQgdGhlIHZlbG9jaXR5IGFwcGxpZWQgdG8gdGhlIGJhbGwgKHRyYWNlIDEpLiBUaGUgeSBkaW1lbnRpb24gb2YgdGhlIHRocm93IGlzIGVzc2VudGlhbGx5IGlnbm9yZWQgKGluIHJlYWxpdHkgdGhlcmUgaXMgYSBzbGlnaHQgdGlsdCBhZGRlZCB0byBhY2NvdW50IGZvciB0aGUgdGlsdCBvZiB0aGUgc3VyZmFjZSkuDQpgYGB7cn0NCnRlc3RfcHB0IDwtIDMNCg0KdGVzdF9kZiA8LSBvbW5pYnVzX2RmICU+JSBmaWx0ZXIocHBpZCA9PSB0ZXN0X3BwdCkNCg0KDQp0cmlhbCA8LSAyNTANCnRyaWFsX2RmIDwtIGZpbHRlcih0ZXN0X2RmLCB0cmlhbF9udW0gPT0gdHJpYWwpDQoNCnggPC0gdHJpYWxfZGYkZmxpY2tfdmVsb2NpdHlfeA0KeSA8LSB0cmlhbF9kZiRmbGlja192ZWxvY2l0eV95DQp6IDwtIHRyaWFsX2RmJGZsaWNrX3ZlbG9jaXR5X3oNCg0KeDIgPC0gdHJpYWxfZGYkZmxpY2tfZGlyZWN0aW9uX3ggKiAtMQ0KeTIgPC0gdHJpYWxfZGYkZmxpY2tfZGlyZWN0aW9uX3kgKiAtMQ0KejIgPC0gdHJpYWxfZGYkZmxpY2tfZGlyZWN0aW9uX3ogKiAtMQ0KDQojIHBsb3QgYm90aA0KcCA8LSBwbG90X2x5KHggPSBjKDAsIHgpLCB5ID0gYygwLCB5KSwgeiA9IGMoMCwgeiksIHR5cGUgPSAic2NhdHRlcjNkIiwgbW9kZSA9ICJsaW5lcyIpICU+JQ0KICBhZGRfdHJhY2UoeCA9IGMoMCwgeDIpLCB5ID0gYygwLCB5MiksIHogPSBjKDAsIHoyKSwgdHlwZSA9ICJzY2F0dGVyM2QiLCBtb2RlID0gImxpbmVzIikgJT4lDQogIGxheW91dChzY2VuZSA9IGxpc3QoDQogICAgeGF4aXMgPSBsaXN0KHRpdGxlID0gIngiLCByYW5nZSA9IGMoLTIsIDIpKSwNCiAgICB5YXhpcyA9IGxpc3QodGl0bGUgPSAieSIsIHJhbmdlID0gYygtMSwgMykpLA0KICAgIHpheGlzID0gbGlzdCh0aXRsZSA9ICJ6IiwgcmFuZ2UgPSBjKC0xLCAzKSkNCiAgKSkNCg0KIyBSZW5kZXIgdGhlIHBsb3QNCnANCmBgYA0Kbm90ZTogdGhpcyBpcyBhIHJvdGF0ZWQgdHJpYWwNCg0KIyMjIERpc3RyaWJ1dGlvbiBvZiBlcnJvcnMNCmBgYHtyfQ0KIyBwbG90IGRpc3RyaWJ1dGlvbiBvZiBlcnJvcl9zaXplDQpwIDwtIGdncGxvdChvbW5pYnVzX2RmLCBhZXMoDQogIHggPSBlcnJvcl9zaXplLA0KICBmaWxsID0gdHlwZQ0KKSkgKw0KICBnZW9tX2hpc3RvZ3JhbShiaW53aWR0aCA9IC41LCBhbHBoYSA9IC42KSArDQogIHRoZW1lX21pbmltYWwoKSArDQogIHRoZW1lKHRleHQgPSBlbGVtZW50X3RleHQoc2l6ZSA9IDExKSkgKw0KICBzY2FsZV9maWxsX21hbnVhbCh2YWx1ZXMgPSBjKCIjZjk5ODJjIiwgIiNkNDAwMDAiKSkgKw0KICBsYWJzKHggPSAiRXJyb3IgU2l6ZSAoY20pIiwgeSA9ICJDb3VudCIpDQoNCnANCmBgYA0KDQojIyMgRGlzdHJpYnV0aW9uIG9mIHRocm93IGFuZ2xlcw0KYGBge3J9DQojIHBsb3QgZGlzdHJpYnV0aW9uIG9mIGVycm9yX3NpemUNCnAgPC0gZ2dwbG90KG9tbmlidXNfZGYsIGFlcygNCiAgeCA9IHRocm93X2RldmlhdGlvbiwNCiAgZmlsbCA9IHR5cGUNCikpICsNCiAgZ2VvbV9oaXN0b2dyYW0oYmlud2lkdGggPSAxLCBhbHBoYSA9IC42KSArDQogIHRoZW1lX21pbmltYWwoKSArDQogIHRoZW1lKHRleHQgPSBlbGVtZW50X3RleHQoc2l6ZSA9IDExKSkgKw0KICBzY2FsZV9maWxsX21hbnVhbCh2YWx1ZXMgPSBjKCIjZjk5ODJjIiwgIiNkNDAwMDAiKSkgKw0KICBsYWJzKA0KICAgIHggPSAiVGhyb3cgQW5nbGUgKMKwKSIsIHkgPSAiQ291bnQiDQogICkgKyAjIGRhc2hlZCBsaW5lcyBhdCAwLCAtMTUsIC0zMA0KICBnZW9tX3ZsaW5lKA0KICAgIHhpbnRlcmNlcHQgPSBjKDAsIC0xNSwgLTMwKSwgbGluZXdpZHRoID0gMC40LA0KICAgIGNvbG91ciA9ICIjQ0NDQ0NDIiwgbGluZXR5cGUgPSAiZGFzaGVkIg0KICApICsgIyB0aWNrcyBvZiAxNSBkZWdyZWVzDQogIHNjYWxlX3hfY29udGludW91cygNCiAgICBicmVha3MgPSBjKC0zMCwgMCwgMzAsIC02MCwgLTkwKQ0KICApDQoNCnANCmBgYA0KDQojIE9yaWdpbmFsIEV4cGVyaW1lbnRzDQoNCiMjIyBQbG90IEFOR1VMQVIgREVWSUFUSU9OUyAoaGFuZCBhbmdsZXMpDQpOb3RlOiBCbHVlcyA9IEFjY2VsZXJhdGlvbiBQZXJ0dXJiYXRpb25zDQpgYGB7ciwgZmlnLndpZHRoPTEwLCBmaWcuaGVpZ2h0PTYsIG91dC53aWR0aD0iMTAwJSJ9DQojIHJlc3Qgb2YgdGhlIGV4cHMNCmRhdGFfcGVyX2dyb3VwIDwtIG9tbmlidXNfZGYgJT4lDQogIGZpbHRlcihleHBfbGFiZWwgPT0gIm9yaWdpbmFsX2V4cHMiIHwgZXhwX2xhYmVsID09ICJjdXJ2ZWRfcGF0aCIpICU+JQ0KICBncm91cF9ieShleHBlcmltZW50LCB0ZXN0X3R5cGUsIHRyaWFsX251bSkgJT4lDQogIHN1bW1hcmlzZSgNCiAgICBtZWFuX2RldmlhdGlvbiA9IG1lYW4odGhyb3dfZGV2aWF0aW9uKSwNCiAgICBjaV9kZXZpYXRpb24gPSB2ZWN0b3JfY29uZmludCh0aHJvd19kZXZpYXRpb24pLA0KICAgIC5ncm91cHMgPSAiZHJvcCINCiAgKQ0KDQojIHNldCB1cCBwbG90DQpwIDwtIGRhdGFfcGVyX2dyb3VwICU+JQ0KICBnZ3Bsb3QoDQogICAgYWVzKA0KICAgICAgeCA9IHRyaWFsX251bSwgeSA9IG1lYW5fZGV2aWF0aW9uLCBjb2xvdXIgPSBleHBlcmltZW50DQogICAgKQ0KICApICsNCiAgdGhlbWVfY2xhc3NpYygpICsNCiAgIyB0aGVtZShsZWdlbmQucG9zaXRpb24gPSAibm9uZSIpICsNCiAgbGFicygNCiAgICB4ID0gIlRyaWFsIE51bWJlciIsDQogICAgeSA9ICJUaHJvdyBBbmdsZSAowrApIg0KICApDQoNCiMgYWRkIGhvcml6b250YWwgbGluZXMNCnAgPC0gcCArDQogIGdlb21faGxpbmUoDQogICAgeWludGVyY2VwdCA9IGMoMCwgLTMwKSwgbGluZXdpZHRoID0gMC40LA0KICAgIGNvbG91ciA9ICIjQ0NDQ0NDIiwgbGluZXR5cGUgPSAic29saWQiDQogICkgKw0KICBnZW9tX2hsaW5lKA0KICAgIHlpbnRlcmNlcHQgPSBjKC0xNSksIGxpbmV3aWR0aCA9IDAuNCwNCiAgICBjb2xvdXIgPSAiI0NDQ0NDQyIsIGxpbmV0eXBlID0gImRhc2hlZCINCiAgKQ0KDQojIHAgPC0gcCArDQojICAgc2NhbGVfeV9jb250aW51b3VzKA0KIyAgICAgbGltaXRzID0gYygtMTAsIDM1KSwNCiMgICAgIGJyZWFrcyA9IGMoMCwgMTUsIDMwKSwNCiMgICAgIGxhYmVscyA9IGMoMCwgMTUsIDMwKQ0KIyAgICkgKw0KIyAgIHNjYWxlX3hfY29udGludW91cygNCiMgICAgIGxpbWl0cyA9IGMoMCwgMTgwKSwNCiMgICAgIGJyZWFrcyA9IGMoMCwgNjAsIDEyMCwgMTgwKSwNCiMgICAgIGxhYmVscyA9IGMoMCwgNjAsIDEyMCwgMTgwKQ0KIyAgICkNCg0KIyBzZXQgZm9udCBzaXplIHRvIDExDQpwIDwtIHAgKw0KICB0aGVtZSh0ZXh0ID0gZWxlbWVudF90ZXh0KHNpemUgPSAxMSkpDQoNCiMgYWRkIGNvbmZpZGVuY2UgaW50ZXJ2YWxzIGFuZCBkYXRhIHBvaW50cw0KcCA8LSBwICsgZ2VvbV9yaWJib24oDQogIGFlcygNCiAgICB5bWluID0gbWVhbl9kZXZpYXRpb24gLSBjaV9kZXZpYXRpb24sDQogICAgeW1heCA9IG1lYW5fZGV2aWF0aW9uICsgY2lfZGV2aWF0aW9uLA0KICAgIGZpbGwgPSBleHBlcmltZW50DQogICksDQogIGNvbG91ciA9IE5BLCBhbHBoYSA9IDAuMw0KKSArIGdlb21fbGluZSgpDQoNCiMgc2V0IGNvbG91ciBwYWxldHRlDQpwIDwtIHAgKw0KICBzY2FsZV9jb2xvdXJfbWFudWFsKHZhbHVlcyA9IHBhbGxldGVfbGlzdCkgKw0KICBzY2FsZV9maWxsX21hbnVhbCh2YWx1ZXMgPSBwYWxsZXRlX2xpc3QpDQoNCg0KIyAjIHNhdmUNCiMgaWYgKHNhdmVfcGxvdHMpIHsNCiMgICBnZ3NhdmUoDQojICAgcCwNCiMgICBmaWxlbmFtZSA9ICIuLi9wbG90cy9wYXBlcl9maWdzL3NyXzMwX3RyYWluaW5nLnBkZiIsIGRldmljZSA9ICJwZGYiLA0KIyAgIGhlaWdodCA9IDQsIHdpZHRoID0gNg0KIyAgICkNCiMgICB9DQoNCmdncGxvdGx5KHApDQoNCiMgcA0KYGBgDQojIyMgVHJpYWwgc2V0cyBvZiBpbnRlcmVzdCBvbmx5DQpOb3RlOiBCbHVlcyA9IEFjY2VsZXJhdGlvbiBQZXJ0dXJiYXRpb25zDQpgYGB7ciwgZmlnLndpZHRoPTEwLCBmaWcuaGVpZ2h0PTYsIG91dC53aWR0aD0iMTAwJSJ9DQojIGZpbHRlciBvdXQganVzdCB0aGUgdHJpYWxzIG9mIGludGVyZXN0DQpkYXRhX3Blcl9ncm91cCA8LSBkYXRhX3Blcl9ncm91cCAlPiUNCiAgZmlsdGVyKA0KICAgIHRlc3RfdHlwZSAhPSAib3RoZXIiDQogICkNCiMgYWRkIGEgZHVtbXkgY29sdW1uIHdpdGggcmVwZWF0aW5nIHNlcXVlbmNlDQojIE5PVEU6IHRoaXMgY2FuJ3QgYmUgY29tYmluZWQgd2l0aCBhYm92ZSBzaW5jZSB3ZSBhcmUgdXNpbmcgbnJvdw0KZGF0YV9wZXJfZ3JvdXAgPC0gZGF0YV9wZXJfZ3JvdXAgJT4lDQogIG11dGF0ZShkdW1teV94ID0gcmVwKDE6KG5yb3coZGF0YV9wZXJfZ3JvdXApIC8gTlVNX0VYUFMpLA0KICAgIGxlbmd0aC5vdXQgPSBucm93KGRhdGFfcGVyX2dyb3VwKQ0KICApKQ0KDQojIHNldCB1cCBwbG90DQpwIDwtIGRhdGFfcGVyX2dyb3VwICU+JQ0KICBnZ3Bsb3QoDQogICAgYWVzKA0KICAgICAgeCA9IGR1bW15X3gsIHkgPSBtZWFuX2RldmlhdGlvbiwgY29sb3VyID0gZXhwZXJpbWVudA0KICAgICkNCiAgKSArDQogIHRoZW1lX2NsYXNzaWMoKSArDQogICMgdGhlbWUobGVnZW5kLnBvc2l0aW9uID0gIm5vbmUiKSArDQogIGxhYnMoDQogICAgeCA9ICJUcmlhbCBOdW1iZXIiLA0KICAgIHkgPSAiVGhyb3cgQW5nbGUgKMKwKSINCiAgKQ0KDQojIGFkZCBob3Jpem9udGFsIGxpbmVzDQpwIDwtIHAgKw0KICBnZW9tX2hsaW5lKA0KICAgIHlpbnRlcmNlcHQgPSBjKDAsIC0zMCksIGxpbmV3aWR0aCA9IDAuNCwNCiAgICBjb2xvdXIgPSAiI0NDQ0NDQyIsIGxpbmV0eXBlID0gInNvbGlkIg0KICApICsNCiAgZ2VvbV9obGluZSgNCiAgICB5aW50ZXJjZXB0ID0gYygtMTUpLCBsaW5ld2lkdGggPSAwLjQsDQogICAgY29sb3VyID0gIiNDQ0NDQ0MiLCBsaW5ldHlwZSA9ICJkYXNoZWQiDQogICkNCg0KIyBhZGQgY29uZmlkZW5jZSBpbnRlcnZhbHMgYW5kIGRhdGEgcG9pbnRzDQpmb3IgKHVuaXF1ZV90ZXN0X3R5cGUgaW4gdW5pcXVlKGRhdGFfcGVyX2dyb3VwJHRlc3RfdHlwZSkpIHsNCiAgIyBnZXQgdGhlIGRhdGEgZm9yIHRoaXMgYmxvY2sNCiAgdG9fcGxvdF9kYXRhIDwtIGZpbHRlcihkYXRhX3Blcl9ncm91cCwgdGVzdF90eXBlID09IHVuaXF1ZV90ZXN0X3R5cGUpDQoNCiAgcCA8LSBwICsgZ2VvbV9yaWJib24oDQogICAgZGF0YSA9IHRvX3Bsb3RfZGF0YSwNCiAgICBhZXMoDQogICAgICB5bWluID0gbWVhbl9kZXZpYXRpb24gLSBjaV9kZXZpYXRpb24sDQogICAgICB5bWF4ID0gbWVhbl9kZXZpYXRpb24gKyBjaV9kZXZpYXRpb24sDQogICAgICBmaWxsID0gZXhwZXJpbWVudA0KICAgICksIGNvbG91ciA9IE5BLCBhbHBoYSA9IDAuMw0KICApICsgZ2VvbV9saW5lKA0KICAgIGRhdGEgPSB0b19wbG90X2RhdGENCiAgKQ0KfQ0KDQojIHNldCBjb2xvdXIgcGFsZXR0ZQ0KcCA8LSBwICsNCiAgc2NhbGVfY29sb3VyX21hbnVhbCh2YWx1ZXMgPSBwYWxsZXRlX2xpc3QpICsNCiAgc2NhbGVfZmlsbF9tYW51YWwodmFsdWVzID0gcGFsbGV0ZV9saXN0KQ0KDQpnZ3Bsb3RseShwKQ0KYGBgDQojIyMgUGxvdCBOT1JNQUxJWkVEIEFOR1VMQVIgREVWSUFUSU9OUyAoaGFuZCBhbmdsZXMpDQpOb3RlOiBCbHVlcyA9IEFjY2VsZXJhdGlvbiBQZXJ0dXJiYXRpb25zDQpgYGB7ciwgZmlnLndpZHRoPTEwLCBmaWcuaGVpZ2h0PTYsIG91dC53aWR0aD0iMTAwJSJ9DQojIHJlc3Qgb2YgdGhlIGV4cHMNCmRhdGFfcGVyX2dyb3VwIDwtIG9tbmlidXNfZGYgJT4lDQogIGZpbHRlcihleHBfbGFiZWwgPT0gIm9yaWdpbmFsX2V4cHMiIHwgZXhwX2xhYmVsID09ICJjdXJ2ZWRfcGF0aCIpICU+JQ0KICBncm91cF9ieShleHBlcmltZW50LCB0ZXN0X3R5cGUsIHRyaWFsX251bSkgJT4lDQogIHN1bW1hcmlzZSgNCiAgICBtZWFuX2RldmlhdGlvbiA9IG1lYW4obm9ybV90aHJvd19kZXZpYXRpb24pLA0KICAgIGNpX2RldmlhdGlvbiA9IHZlY3Rvcl9jb25maW50KG5vcm1fdGhyb3dfZGV2aWF0aW9uKSwNCiAgICAuZ3JvdXBzID0gImRyb3AiDQogICkNCg0KIyBzZXQgdXAgcGxvdA0KcCA8LSBkYXRhX3Blcl9ncm91cCAlPiUNCiAgZ2dwbG90KA0KICAgIGFlcygNCiAgICAgIHggPSB0cmlhbF9udW0sIHkgPSBtZWFuX2RldmlhdGlvbiwgY29sb3VyID0gZXhwZXJpbWVudA0KICAgICkNCiAgKSArDQogIHRoZW1lX2NsYXNzaWMoKSArDQogICMgdGhlbWUobGVnZW5kLnBvc2l0aW9uID0gIm5vbmUiKSArDQogIGxhYnMoDQogICAgeCA9ICJUcmlhbCBOdW1iZXIiLA0KICAgIHkgPSAiTm9ybWFsaXplZCBUaHJvdyBBbmdsZSINCiAgKQ0KDQojIGFkZCBob3Jpem9udGFsIGxpbmVzDQpwIDwtIHAgKw0KICBnZW9tX2hsaW5lKA0KICAgIHlpbnRlcmNlcHQgPSBjKDAsIDEsIDIpLCBsaW5ld2lkdGggPSAwLjQsDQogICAgY29sb3VyID0gIiNDQ0NDQ0MiLCBsaW5ldHlwZSA9ICJzb2xpZCINCiAgKSArDQogIGdlb21faGxpbmUoDQogICAgeWludGVyY2VwdCA9IGMoMC41LCAxLjUpLCBsaW5ld2lkdGggPSAwLjQsDQogICAgY29sb3VyID0gIiNDQ0NDQ0MiLCBsaW5ldHlwZSA9ICJkYXNoZWQiDQogICkNCg0KIyBwIDwtIHAgKw0KIyAgIHNjYWxlX3lfY29udGludW91cygNCiMgICAgIGxpbWl0cyA9IGMoLTEwLCAzNSksDQojICAgICBicmVha3MgPSBjKDAsIDE1LCAzMCksDQojICAgICBsYWJlbHMgPSBjKDAsIDE1LCAzMCkNCiMgICApICsNCiMgICBzY2FsZV94X2NvbnRpbnVvdXMoDQojICAgICBsaW1pdHMgPSBjKDAsIDE4MCksDQojICAgICBicmVha3MgPSBjKDAsIDYwLCAxMjAsIDE4MCksDQojICAgICBsYWJlbHMgPSBjKDAsIDYwLCAxMjAsIDE4MCkNCiMgICApDQoNCiMgc2V0IGZvbnQgc2l6ZSB0byAxMQ0KcCA8LSBwICsNCiAgdGhlbWUodGV4dCA9IGVsZW1lbnRfdGV4dChzaXplID0gMTEpKQ0KDQojIGFkZCBjb25maWRlbmNlIGludGVydmFscyBhbmQgZGF0YSBwb2ludHMNCnAgPC0gcCArIGdlb21fcmliYm9uKA0KICBhZXMoDQogICAgeW1pbiA9IG1lYW5fZGV2aWF0aW9uIC0gY2lfZGV2aWF0aW9uLA0KICAgIHltYXggPSBtZWFuX2RldmlhdGlvbiArIGNpX2RldmlhdGlvbiwNCiAgICBmaWxsID0gZXhwZXJpbWVudA0KICApLA0KICBjb2xvdXIgPSBOQSwgYWxwaGEgPSAwLjMNCikgKyBnZW9tX2xpbmUoKQ0KDQojIHNldCBjb2xvdXIgcGFsZXR0ZQ0KcCA8LSBwICsNCiAgc2NhbGVfY29sb3VyX21hbnVhbCh2YWx1ZXMgPSBwYWxsZXRlX2xpc3QpICsNCiAgc2NhbGVfZmlsbF9tYW51YWwodmFsdWVzID0gcGFsbGV0ZV9saXN0KQ0KDQoNCiMgIyBzYXZlDQojIGlmIChzYXZlX3Bsb3RzKSB7DQojICAgZ2dzYXZlKA0KIyAgIHAsDQojICAgZmlsZW5hbWUgPSAiLi4vcGxvdHMvcGFwZXJfZmlncy9zcl8zMF90cmFpbmluZy5wZGYiLCBkZXZpY2UgPSAicGRmIiwNCiMgICBoZWlnaHQgPSA0LCB3aWR0aCA9IDYNCiMgICApDQojICAgfQ0KDQpnZ3Bsb3RseShwKQ0KDQojIHANCmBgYA0KIyMjIFRyaWFsIHNldHMgb2YgaW50ZXJlc3Qgb25seQ0KTm90ZTogQmx1ZXMgPSBBY2NlbGVyYXRpb24gUGVydHVyYmF0aW9ucw0KYGBge3IsIGZpZy53aWR0aD0xMCwgZmlnLmhlaWdodD02LCBvdXQud2lkdGg9IjEwMCUifQ0KIyBmaWx0ZXIgb3V0IGp1c3QgdGhlIHRyaWFscyBvZiBpbnRlcmVzdA0KZGF0YV9wZXJfZ3JvdXAgPC0gZGF0YV9wZXJfZ3JvdXAgJT4lDQogIGZpbHRlcigNCiAgICB0ZXN0X3R5cGUgIT0gIm90aGVyIg0KICApDQojIGFkZCBhIGR1bW15IGNvbHVtbiB3aXRoIHJlcGVhdGluZyBzZXF1ZW5jZQ0KIyBOT1RFOiB0aGlzIGNhbid0IGJlIGNvbWJpbmVkIHdpdGggYWJvdmUgc2luY2Ugd2UgYXJlIHVzaW5nIG5yb3cNCmRhdGFfcGVyX2dyb3VwIDwtIGRhdGFfcGVyX2dyb3VwICU+JQ0KICBtdXRhdGUoZHVtbXlfeCA9IHJlcCgxOihucm93KGRhdGFfcGVyX2dyb3VwKSAvIE5VTV9FWFBTKSwNCiAgICBsZW5ndGgub3V0ID0gbnJvdyhkYXRhX3Blcl9ncm91cCkNCiAgKSkNCg0KIyBzZXQgdXAgcGxvdA0KcCA8LSBkYXRhX3Blcl9ncm91cCAlPiUNCiAgZ2dwbG90KA0KICAgIGFlcygNCiAgICAgIHggPSBkdW1teV94LCB5ID0gbWVhbl9kZXZpYXRpb24sIGNvbG91ciA9IGV4cGVyaW1lbnQNCiAgICApDQogICkgKw0KICB0aGVtZV9jbGFzc2ljKCkgKw0KICAjIHRoZW1lKGxlZ2VuZC5wb3NpdGlvbiA9ICJub25lIikgKw0KICBsYWJzKA0KICAgIHggPSAiVHJpYWwgTnVtYmVyIiwNCiAgICB5ID0gIk5vcm1hbGl6ZWQgVGhyb3cgQW5nbGUiDQogICkNCg0KIyBhZGQgaG9yaXpvbnRhbCBsaW5lcw0KcCA8LSBwICsNCiAgZ2VvbV9obGluZSgNCiAgICB5aW50ZXJjZXB0ID0gYygwLCAxLCAyKSwgbGluZXdpZHRoID0gMC40LA0KICAgIGNvbG91ciA9ICIjQ0NDQ0NDIiwgbGluZXR5cGUgPSAic29saWQiDQogICkgKw0KICBnZW9tX2hsaW5lKA0KICAgIHlpbnRlcmNlcHQgPSBjKDAuNSwgMS41KSwgbGluZXdpZHRoID0gMC40LA0KICAgIGNvbG91ciA9ICIjQ0NDQ0NDIiwgbGluZXR5cGUgPSAiZGFzaGVkIg0KICApDQoNCiMgYWRkIGNvbmZpZGVuY2UgaW50ZXJ2YWxzIGFuZCBkYXRhIHBvaW50cw0KZm9yICh1bmlxdWVfdGVzdF90eXBlIGluIHVuaXF1ZShkYXRhX3Blcl9ncm91cCR0ZXN0X3R5cGUpKSB7DQogICMgZ2V0IHRoZSBkYXRhIGZvciB0aGlzIGJsb2NrDQogIHRvX3Bsb3RfZGF0YSA8LSBmaWx0ZXIoZGF0YV9wZXJfZ3JvdXAsIHRlc3RfdHlwZSA9PSB1bmlxdWVfdGVzdF90eXBlKQ0KDQogIHAgPC0gcCArIGdlb21fcmliYm9uKA0KICAgIGRhdGEgPSB0b19wbG90X2RhdGEsDQogICAgYWVzKA0KICAgICAgeW1pbiA9IG1lYW5fZGV2aWF0aW9uIC0gY2lfZGV2aWF0aW9uLA0KICAgICAgeW1heCA9IG1lYW5fZGV2aWF0aW9uICsgY2lfZGV2aWF0aW9uLA0KICAgICAgZmlsbCA9IGV4cGVyaW1lbnQNCiAgICApLCBjb2xvdXIgPSBOQSwgYWxwaGEgPSAwLjMNCiAgKSArIGdlb21fbGluZSgNCiAgICBkYXRhID0gdG9fcGxvdF9kYXRhDQogICkNCn0NCg0KIyBzZXQgY29sb3VyIHBhbGV0dGUNCnAgPC0gcCArDQogIHNjYWxlX2NvbG91cl9tYW51YWwodmFsdWVzID0gcGFsbGV0ZV9saXN0KSArDQogIHNjYWxlX2ZpbGxfbWFudWFsKHZhbHVlcyA9IHBhbGxldGVfbGlzdCkNCg0KZ2dwbG90bHkocCkNCmBgYA0KDQojIyBMZWFybmluZyBSYXRlcywgU3RhcnQgUG9pbnRzIGFuZCBBc3ltcHRvdGVzDQojIyMgTGVhcm5pbmcgUmF0ZXMNCmBgYHtyLCBmaWcud2lkdGg9MTAsIGZpZy5oZWlnaHQ9Niwgb3V0LndpZHRoPSIxMDAlIn0NCmRhdGFfcHB0IDwtIGJsb2NrX2luaXRfbGVhcm5pbmdfcmF0ZXMgJT4lDQogIGZpbHRlcihleHBlcmltZW50ICE9ICJhX2JhbGxfcm9sbF9hbmltYXRlX3N1cmZhY2UiKQ0KDQpkYXRhX2dyb3VwIDwtIGRhdGFfcHB0ICU+JQ0KICBncm91cF9ieShleHBlcmltZW50LCB0ZXN0X3R5cGUpICU+JQ0KICBzdW1tYXJpc2UoDQogICAgbWVhbl9sZWFybmluZ19yYXRlID0gbWVhbihleHBfZml0X2xhbWJkYSksDQogICAgY2lfbGVhcm5pbmdfcmF0ZSA9IHZlY3Rvcl9jb25maW50KGV4cF9maXRfbGFtYmRhKSwNCiAgICBtZWFuX2hpZ2ggPSBtZWFuKGV4cF9maXRfTjApLA0KICAgIGNpX2hpZ2ggPSB2ZWN0b3JfY29uZmludChleHBfZml0X04wKSwNCiAgICBtZWFuX2xvdyA9IG1lYW4oZXhwX2ZpdF9kaXNwbGFjZSksDQogICAgY2lfbG93ID0gdmVjdG9yX2NvbmZpbnQoZXhwX2ZpdF9kaXNwbGFjZSksDQogICAgLmdyb3VwcyA9ICJkcm9wIg0KICApDQoNCnAgPC0gZGF0YV9ncm91cCAlPiUNCiAgZ2dwbG90KA0KICAgIGFlcyh4ID0gZXhwZXJpbWVudCwgeSA9IG1lYW5fbGVhcm5pbmdfcmF0ZSwgY29sb3VyID0gZXhwZXJpbWVudCkNCiAgKSArDQogIHRoZW1lX2NsYXNzaWMoKSArDQogIGxhYnMoDQogICAgeCA9IE5VTEwsDQogICAgeSA9ICJMZWFybmluZyBSYXRlIg0KICApICsNCiAgZmFjZXRfd3JhcCh+dGVzdF90eXBlKQ0KDQojIHJlbW92ZSBhbGwgeCBheGlzIGxhYmVscw0KcCA8LSBwICsgdGhlbWUoYXhpcy50ZXh0LnggPSBlbGVtZW50X2JsYW5rKCkpDQoNCiMgZm9yIHRoZSBjb2xvdXIgbGVnZW5kLCBvbmx5IHNob3cgdGhlIGZpcnN0IDcgTm90ZSB0aGlzIGRvZXNuJ3Qgd29yayBmb3IgdGhlIHBsb3RseSBwbG90DQpwIDwtIHAgKyBndWlkZXMoY29sb3VyID0gZ3VpZGVfbGVnZW5kKG92ZXJyaWRlLmFlcyA9IGxpc3QoYWxwaGEgPSAxKSkpDQoNCiMgYWRkIGRhdGEgcG9pbnRzDQpwIDwtIHAgKw0KICBnZW9tX2JlZXN3YXJtKA0KICAgIGRhdGEgPSBkYXRhX3BwdCwNCiAgICBhZXMoDQogICAgICB5ID0gZXhwX2ZpdF9sYW1iZGENCiAgICApLA0KICAgIGFscGhhID0gMC4xLA0KICAgIHNpemUgPSAxDQogICkgKw0KICBnZW9tX2xpbmVyYW5nZShhZXMoDQogICAgeW1pbiA9IG1lYW5fbGVhcm5pbmdfcmF0ZSAtIGNpX2xlYXJuaW5nX3JhdGUsDQogICAgeW1heCA9IG1lYW5fbGVhcm5pbmdfcmF0ZSArIGNpX2xlYXJuaW5nX3JhdGUNCiAgKSwgYWxwaGEgPSAwLjUsIGx3ZCA9IDIpICsNCiAgZ2VvbV9wb2ludCgpDQoNCiMgc2V0IGNvbG91ciBwYWxldHRlDQpwIDwtIHAgKw0KICBzY2FsZV9jb2xvdXJfbWFudWFsKHZhbHVlcyA9IHBhbGxldGVfbGlzdCkgKw0KICBzY2FsZV9maWxsX21hbnVhbCh2YWx1ZXMgPSBwYWxsZXRlX2xpc3QpDQoNCmdncGxvdGx5KHApDQojIHANCmBgYA0KDQojIyMgSGlnaCBQb2ludHMgKEFzeW1wdG90ZXMgb3IgU3RhcnRzKQ0KYGBge3IsIGZpZy53aWR0aD0xMCwgZmlnLmhlaWdodD02LCBvdXQud2lkdGg9IjEwMCUifQ0KcCA8LSBkYXRhX2dyb3VwICU+JQ0KICBnZ3Bsb3QoDQogICAgYWVzKHggPSBleHBlcmltZW50LCB5ID0gbWVhbl9oaWdoLCBjb2xvdXIgPSBleHBlcmltZW50KQ0KICApICsNCiAgdGhlbWVfY2xhc3NpYygpICsNCiAgbGFicygNCiAgICB4ID0gTlVMTCwNCiAgICB5ID0gIkhpZ2ggUG9pbnQiDQogICkgKw0KICBmYWNldF93cmFwKH50ZXN0X3R5cGUpDQoNCiMgcmVtb3ZlIGFsbCB4IGF4aXMgbGFiZWxzDQpwIDwtIHAgKyB0aGVtZShheGlzLnRleHQueCA9IGVsZW1lbnRfYmxhbmsoKSkNCg0KIyBmb3IgdGhlIGNvbG91ciBsZWdlbmQsIG9ubHkgc2hvdyB0aGUgZmlyc3QgNyBOb3RlIHRoaXMgZG9lc24ndCB3b3JrIGZvciB0aGUgcGxvdGx5IHBsb3QNCnAgPC0gcCArIGd1aWRlcyhjb2xvdXIgPSBndWlkZV9sZWdlbmQob3ZlcnJpZGUuYWVzID0gbGlzdChhbHBoYSA9IDEpKSkNCg0KIyBhZGQgZGF0YSBwb2ludHMNCnAgPC0gcCArDQogIGdlb21fYmVlc3dhcm0oDQogICAgZGF0YSA9IGRhdGFfcHB0LA0KICAgIGFlcygNCiAgICAgIHkgPSBleHBfZml0X04wDQogICAgKSwNCiAgICBhbHBoYSA9IDAuMSwNCiAgICBzaXplID0gMQ0KICApICsNCiAgZ2VvbV9saW5lcmFuZ2UoYWVzKA0KICAgIHltaW4gPSBtZWFuX2hpZ2ggLSBjaV9oaWdoLA0KICAgIHltYXggPSBtZWFuX2hpZ2ggKyBjaV9oaWdoDQogICksIGFscGhhID0gMC41LCBsd2QgPSAyKSArDQogIGdlb21fcG9pbnQoKQ0KDQojIHNldCBjb2xvdXIgcGFsZXR0ZQ0KcCA8LSBwICsNCiAgc2NhbGVfY29sb3VyX21hbnVhbCh2YWx1ZXMgPSBwYWxsZXRlX2xpc3QpICsNCiAgc2NhbGVfZmlsbF9tYW51YWwodmFsdWVzID0gcGFsbGV0ZV9saXN0KQ0KDQpnZ3Bsb3RseShwKQ0KIyBwDQpgYGANCiMjIyBMb3cgUG9pbnRzIChBc3ltcHRvdGVzIG9yIFN0YXJ0cykNCmBgYHtyLCBmaWcud2lkdGg9MTAsIGZpZy5oZWlnaHQ9Niwgb3V0LndpZHRoPSIxMDAlIn0NCnAgPC0gZGF0YV9ncm91cCAlPiUNCiAgZ2dwbG90KA0KICAgIGFlcyh4ID0gZXhwZXJpbWVudCwgeSA9IG1lYW5fbG93LCBjb2xvdXIgPSBleHBlcmltZW50KQ0KICApICsNCiAgdGhlbWVfY2xhc3NpYygpICsNCiAgbGFicygNCiAgICB4ID0gTlVMTCwNCiAgICB5ID0gIkxvdyBQb2ludCINCiAgKSArDQogIGZhY2V0X3dyYXAofnRlc3RfdHlwZSkNCg0KIyByZW1vdmUgYWxsIHggYXhpcyBsYWJlbHMNCnAgPC0gcCArIHRoZW1lKGF4aXMudGV4dC54ID0gZWxlbWVudF9ibGFuaygpKQ0KDQojIGZvciB0aGUgY29sb3VyIGxlZ2VuZCwgb25seSBzaG93IHRoZSBmaXJzdCA3IE5vdGUgdGhpcyBkb2Vzbid0IHdvcmsgZm9yIHRoZSBwbG90bHkgcGxvdA0KcCA8LSBwICsgZ3VpZGVzKGNvbG91ciA9IGd1aWRlX2xlZ2VuZChvdmVycmlkZS5hZXMgPSBsaXN0KGFscGhhID0gMSkpKQ0KDQojIGFkZCBkYXRhIHBvaW50cw0KcCA8LSBwICsNCiAgZ2VvbV9iZWVzd2FybSgNCiAgICBkYXRhID0gZGF0YV9wcHQsDQogICAgYWVzKA0KICAgICAgeSA9IGV4cF9maXRfZGlzcGxhY2UNCiAgICApLA0KICAgIGFscGhhID0gMC4xLA0KICAgIHNpemUgPSAxDQogICkgKw0KICBnZW9tX2xpbmVyYW5nZShhZXMoDQogICAgeW1pbiA9IG1lYW5fbG93IC0gY2lfbG93LA0KICAgIHltYXggPSBtZWFuX2xvdyArIGNpX2xvdw0KICApLCBhbHBoYSA9IDAuNSwgbHdkID0gMikgKw0KICBnZW9tX3BvaW50KCkNCg0KIyBzZXQgY29sb3VyIHBhbGV0dGUNCnAgPC0gcCArDQogIHNjYWxlX2NvbG91cl9tYW51YWwodmFsdWVzID0gcGFsbGV0ZV9saXN0KSArDQogIHNjYWxlX2ZpbGxfbWFudWFsKHZhbHVlcyA9IHBhbGxldGVfbGlzdCkNCg0KZ2dwbG90bHkocCkNCiMgcA0KYGBgDQoNCkZvciB3YXNob3V0OiBUaGUgQ1VFRCBhY2NlbCArIGN1cnZlZCBoYXZlIGEgbG93ZXIgc3RhcnRpbmcgcG9pbnQgKHRoZXJlZm9yZSAtLSBjdWUgd29ya3MpLiBWTVIgZ3JvdXAgaGFzIHNsaWdodGx5IGxvd2VyLg0KV2hlbiBjb21wYXJpbmcgZXZlcnl0aGluZyB3aXRoIGEgaGlnaCBzdGFydGluZyBwb2ludCwgdGhlICBBQ0NFTCBncm91cCBoYXMgYSBtdWNoIEZBU1RFUiBsZWFybmluZyByYXRlLiANCldoZW4gdHJhbnNmZXJyaW5nLCBubyBkaWZmZXJlbmNlIGluIGxlYXJuaW5nIHJhdGVzLg0KDQpgYGB7ciwgZmlnLndpZHRoPTEwLCBmaWcuaGVpZ2h0PTYsIG91dC53aWR0aD0iMTAwJSJ9DQojIENvbXBhcmUgbGVhcm5pbmcgcmF0ZXMgYW5kIGhpZ2ggcG9pbnRzDQpwIDwtIGRhdGFfcHB0ICU+JSANCiAgZ2dwbG90KA0KICAgIGFlcyh4ID0gZXhwX2ZpdF9sYW1iZGEsIHkgPSBleHBfZml0X04wLCBjb2xvdXIgPSBleHBlcmltZW50KQ0KICApICsNCiAgdGhlbWVfY2xhc3NpYygpICsNCiAgbGFicygNCiAgICB4ID0gIkxlYXJuaW5nIFJhdGUiLA0KICAgIHkgPSAiSGlnaCBQb2ludCINCiAgKSArDQogIGZhY2V0X3dyYXAofnRlc3RfdHlwZSkgKw0KICBnZW9tX3BvaW50KCkgKw0KICBzY2FsZV9jb2xvdXJfbWFudWFsKHZhbHVlcyA9IHBhbGxldGVfbGlzdCkgKw0KICBzY2FsZV9maWxsX21hbnVhbCh2YWx1ZXMgPSBwYWxsZXRlX2xpc3QpDQoNCmdncGxvdGx5KHApDQpgYGANCg0KIyMjIENvcnJlbGF0aW9ucyB3aXRoaW4gdGhlIDMgdmFyaWFibGVzIGFib3ZlDQoNCmBgYHtyLCBmaWcud2lkdGg9MTAsIGZpZy5oZWlnaHQ9Niwgb3V0LndpZHRoPSIxMDAlIiwgd2FybmluZz1GQUxTRX0NCg0KIyAzZCBwbG90IG9mIGxlYXJuaW5nIHJhdGVzLCBoaWdoIHBvaW50cywgYW5kIGxvdyBwb2ludHMNCiMgZmlsdGVyIG91dCBqdXN0IHRoZSB3YXNob3V0IHRyaWFscyBmb3Igbm93DQpkYXRhX3BwdCA8LSBibG9ja19pbml0X2xlYXJuaW5nX3JhdGVzICU+JQ0KICBmaWx0ZXIodGVzdF90eXBlID09ICJ3YXNob3V0X2luaXQiKQ0KDQpwIDwtIHBsb3RfbHkodHlwZSA9ICdzY2F0dGVyM2QnLCBtb2RlID0gJ21hcmtlcnMnKQ0KDQpmb3IgKGV4cGVyaW1lbnRfIGluIHVuaXF1ZShkYXRhX3BwdCRleHBlcmltZW50KSkgew0KICAjIGdldCB0aGUgZGF0YSBmb3IgdGhpcyBibG9jaw0KICB0b19wbG90X2RhdGEgPC0gZmlsdGVyKGRhdGFfcHB0LCBleHBlcmltZW50ID09IGV4cGVyaW1lbnRfKQ0KICANCiAgcCA8LSBwICU+JQ0KICAgIGFkZF90cmFjZSgNCiAgICAgIHggPSB0b19wbG90X2RhdGEkZXhwX2ZpdF9sYW1iZGEsIA0KICAgICAgeSA9IHRvX3Bsb3RfZGF0YSRleHBfZml0X04wLCANCiAgICAgIHogPSB0b19wbG90X2RhdGEkZXhwX2ZpdF9kaXNwbGFjZSwgDQogICAgICBjb2xvciA9IEkocGFsbGV0ZV9saXN0W2V4cGVyaW1lbnRfXSksDQogICAgICBuYW1lID0gZXhwZXJpbWVudF8NCiAgICApDQp9DQoNCnANCmBgYA0KDQotIFdlIG1pZ2h0IHdhbnQgdG8gdGVzdCBpZiBkYXRhIHBvaW50cyBjbHVzdGVyIGRpZmZlcmVudGx5IGRlcGVuZGluZyBvbiB0aGUgZXhwZXJpbWVudCBjb25kaXRpb25zLg0KLSBOb3RlOiBJdCBpcyBhIGdvb2QgaWRlYSB0byBkbyBzb21lIEVEQSB0byBjaGVjayBpZiBhbGwgMyBkaW1lbnNpb25zIGFyZSByZXF1aXJlZC4gQ29ycmVsYXRpb24gbWF0cml4Pw0KLSANCg0KYGBge3J9DQojIyMgVEVTVElORw0KIyByZXN0IG9mIHRoZSBleHBzDQpkYXRhXyA8LSBvbW5pYnVzX2RmICU+JQ0KICBmaWx0ZXIoZXhwZXJpbWVudCA9PSAiYWNjZWxfdW5jdWVkIiwgdGVzdF90eXBlID09ICJ3YXNob3V0X2luaXQiKQ0KICANCiMgc2V0IHVwIHBsb3QNCnAgPC0gZGF0YV8gJT4lDQogIGdncGxvdCgNCiAgICBhZXMoDQogICAgICB4ID0gdHJpYWxfbnVtLCB5ID0gbm9ybV90aHJvd19kZXZpYXRpb24sIGNvbG91ciA9IHBwaWQNCiAgICApDQogICkgKw0KICB0aGVtZV9jbGFzc2ljKCkgKw0KICAjIHRoZW1lKGxlZ2VuZC5wb3NpdGlvbiA9ICJub25lIikgKw0KICBsYWJzKA0KICAgIHggPSAiVHJpYWwgTnVtYmVyIiwNCiAgICB5ID0gIk5vcm1hbGl6ZWQgVGhyb3cgQW5nbGUiDQogICkNCg0KIyBhZGQgaG9yaXpvbnRhbCBsaW5lcw0KcCA8LSBwICsNCiAgZ2VvbV9obGluZSgNCiAgICB5aW50ZXJjZXB0ID0gYygwLCAxLCAyKSwgbGluZXdpZHRoID0gMC40LA0KICAgIGNvbG91ciA9ICIjQ0NDQ0NDIiwgbGluZXR5cGUgPSAic29saWQiDQogICkgKw0KICBnZW9tX2hsaW5lKA0KICAgIHlpbnRlcmNlcHQgPSBjKDAuNSwgMS41KSwgbGluZXdpZHRoID0gMC40LA0KICAgIGNvbG91ciA9ICIjQ0NDQ0NDIiwgbGluZXR5cGUgPSAiZGFzaGVkIg0KICApDQoNCiMgcCA8LSBwICsNCiMgICBzY2FsZV95X2NvbnRpbnVvdXMoDQojICAgICBsaW1pdHMgPSBjKC0xMCwgMzUpLA0KIyAgICAgYnJlYWtzID0gYygwLCAxNSwgMzApLA0KIyAgICAgbGFiZWxzID0gYygwLCAxNSwgMzApDQojICAgKSArDQojICAgc2NhbGVfeF9jb250aW51b3VzKA0KIyAgICAgbGltaXRzID0gYygwLCAxODApLA0KIyAgICAgYnJlYWtzID0gYygwLCA2MCwgMTIwLCAxODApLA0KIyAgICAgbGFiZWxzID0gYygwLCA2MCwgMTIwLCAxODApDQojICAgKQ0KDQojIHNldCBmb250IHNpemUgdG8gMTENCnAgPC0gcCArDQogIHRoZW1lKHRleHQgPSBlbGVtZW50X3RleHQoc2l6ZSA9IDExKSkNCg0KIyBhZGQgY29uZmlkZW5jZSBpbnRlcnZhbHMgYW5kIGRhdGEgcG9pbnRzDQpwIDwtIHAgKyBnZW9tX2xpbmUoKQ0KDQoNCg0KIyAjIHNhdmUNCiMgaWYgKHNhdmVfcGxvdHMpIHsNCiMgICBnZ3NhdmUoDQojICAgcCwNCiMgICBmaWxlbmFtZSA9ICIuLi9wbG90cy9wYXBlcl9maWdzL3NyXzMwX3RyYWluaW5nLnBkZiIsIGRldmljZSA9ICJwZGYiLA0KIyAgIGhlaWdodCA9IDQsIHdpZHRoID0gNg0KIyAgICkNCiMgICB9DQoNCmdncGxvdGx5KHApDQoNCiMgcA0KYGBgDQoNCiMjIFBsb3QgRVJST1IgU0laRQ0KTm90ZTogQmx1ZXMgPSBBY2NlbGVyYXRpb24gUGVydHVyYmF0aW9ucw0KYGBge3IsIGZpZy53aWR0aD0xMCwgZmlnLmhlaWdodD02LCBvdXQud2lkdGg9IjEwMCUifQ0KIyBvcmlnaW5hbCBleHBlcmltZW50cyBvbmx5DQpkYXRhX3Blcl9ncm91cCA8LSBvbW5pYnVzX2RmICU+JQ0KICBmaWx0ZXIoZXhwX2xhYmVsID09ICJvcmlnaW5hbF9leHBzIiB8IGV4cF9sYWJlbCA9PSAiY3VydmVkX3BhdGgiKSAlPiUNCiAgZ3JvdXBfYnkoZXhwZXJpbWVudCwgdGVzdF90eXBlLCB0cmlhbF9udW0pICU+JQ0KICBzdW1tYXJpc2UoDQogICAgbWVhbl9kZXZpYXRpb24gPSBtZWFuKGVycm9yX3NpemUpLA0KICAgIGNpX2RldmlhdGlvbiA9IHZlY3Rvcl9jb25maW50KGVycm9yX3NpemUpLA0KICAgIC5ncm91cHMgPSAiZHJvcCINCiAgKQ0KDQojIHNldCB1cCBwbG90DQpwIDwtIGRhdGFfcGVyX2dyb3VwICU+JQ0KICBnZ3Bsb3QoDQogICAgYWVzKA0KICAgICAgeCA9IHRyaWFsX251bSwgeSA9IG1lYW5fZGV2aWF0aW9uLCBjb2xvdXIgPSBleHBlcmltZW50DQogICAgKQ0KICApICsNCiAgdGhlbWVfY2xhc3NpYygpICsNCiAgIyB0aGVtZShsZWdlbmQucG9zaXRpb24gPSAibm9uZSIpICsNCiAgbGFicygNCiAgICB4ID0gIlRyaWFsIE51bWJlciIsDQogICAgeSA9ICJBYnNvbHV0ZSBUYXJnZXQgRXJyb3IgKGNtKSINCiAgKQ0KDQojIGFkZCBob3Jpem9udGFsIGxpbmVzDQpwIDwtIHAgKw0KICBnZW9tX2hsaW5lKA0KICAgIHlpbnRlcmNlcHQgPSBjKDAsIDQwKSwgbGluZXdpZHRoID0gMC40LA0KICAgIGNvbG91ciA9ICIjQ0NDQ0NDIiwgbGluZXR5cGUgPSAic29saWQiDQogICkgKw0KICBnZW9tX2hsaW5lKA0KICAgIHlpbnRlcmNlcHQgPSBjKDIwKSwgbGluZXdpZHRoID0gMC40LA0KICAgIGNvbG91ciA9ICIjQ0NDQ0NDIiwgbGluZXR5cGUgPSAiZGFzaGVkIg0KICApDQoNCiMgcCA8LSBwICsNCiMgICBzY2FsZV95X2NvbnRpbnVvdXMoDQojICAgICBsaW1pdHMgPSBjKC0xMCwgMzUpLA0KIyAgICAgYnJlYWtzID0gYygwLCAxNSwgMzApLA0KIyAgICAgbGFiZWxzID0gYygwLCAxNSwgMzApDQojICAgKSArDQojICAgc2NhbGVfeF9jb250aW51b3VzKA0KIyAgICAgbGltaXRzID0gYygwLCAxODApLA0KIyAgICAgYnJlYWtzID0gYygwLCA2MCwgMTIwLCAxODApLA0KIyAgICAgbGFiZWxzID0gYygwLCA2MCwgMTIwLCAxODApDQojICAgKQ0KDQojIHNldCBmb250IHNpemUgdG8gMTENCnAgPC0gcCArDQogIHRoZW1lKHRleHQgPSBlbGVtZW50X3RleHQoc2l6ZSA9IDExKSkNCg0KIyBhZGQgY29uZmlkZW5jZSBpbnRlcnZhbHMgYW5kIGRhdGEgcG9pbnRzDQpwIDwtIHAgKyBnZW9tX3JpYmJvbigNCiAgYWVzKA0KICAgIHltaW4gPSBtZWFuX2RldmlhdGlvbiAtIGNpX2RldmlhdGlvbiwNCiAgICB5bWF4ID0gbWVhbl9kZXZpYXRpb24gKyBjaV9kZXZpYXRpb24sDQogICAgZmlsbCA9IGV4cGVyaW1lbnQNCiAgKSwNCiAgY29sb3VyID0gTkEsIGFscGhhID0gMC4zDQopICsgZ2VvbV9saW5lKCkNCg0KIyBzZXQgY29sb3VyIHBhbGV0dGUNCnAgPC0gcCArDQogIHNjYWxlX2NvbG91cl9tYW51YWwodmFsdWVzID0gcGFsbGV0ZV9saXN0KSArDQogIHNjYWxlX2ZpbGxfbWFudWFsKHZhbHVlcyA9IHBhbGxldGVfbGlzdCkNCg0KDQojICMgc2F2ZQ0KIyBpZiAoc2F2ZV9wbG90cykgew0KIyAgIGdnc2F2ZSgNCiMgICBwLA0KIyAgIGZpbGVuYW1lID0gIi4uL3Bsb3RzL3BhcGVyX2ZpZ3Mvc3JfMzBfdHJhaW5pbmcucGRmIiwgZGV2aWNlID0gInBkZiIsDQojICAgaGVpZ2h0ID0gNCwgd2lkdGggPSA2DQojICAgKQ0KIyAgIH0NCg0KZ2dwbG90bHkocCkNCg0KIyBwDQpgYGANCg0KdmlzaWJsZSB2cyBub24tdmlzaWJsZSB0aWx0IGRvZXNuJ3QgYWZmZWN0IHRoZSAxNS1kZWdyZWUgcm90YXRpb24gY29uZGl0aW9uLiBCdXQgYWZmZWN0cyBhbGwgb3RoZXIgY29uZGl0aW9ucy4gU28gMTUtZGVncmVlIHJvdGF0aW9uDQoNCiMjIyBUcmlhbCBzZXRzIG9mIGludGVyZXN0IG9ubHkNCk5vdGU6IEJsdWVzID0gQWNjZWxlcmF0aW9uIFBlcnR1cmJhdGlvbnMNCmBgYHtyLCBmaWcud2lkdGg9MTAsIGZpZy5oZWlnaHQ9Niwgb3V0LndpZHRoPSIxMDAlIn0NCiMgZmlsdGVyIG91dCBqdXN0IHRoZSB0cmlhbHMgb2YgaW50ZXJlc3QNCmRhdGFfcGVyX2dyb3VwIDwtIGRhdGFfcGVyX2dyb3VwICU+JQ0KICBmaWx0ZXIoDQogICAgdGVzdF90eXBlICE9ICJvdGhlciINCiAgKQ0KIyBhZGQgYSBkdW1teSBjb2x1bW4gd2l0aCByZXBlYXRpbmcgc2VxdWVuY2UNCiMgTk9URTogdGhpcyBjYW4ndCBiZSBjb21iaW5lZCB3aXRoIGFib3ZlIHNpbmNlIHdlIGFyZSB1c2luZyBucm93DQpkYXRhX3Blcl9ncm91cCA8LSBkYXRhX3Blcl9ncm91cCAlPiUNCiAgbXV0YXRlKGR1bW15X3ggPSByZXAoMToobnJvdyhkYXRhX3Blcl9ncm91cCkgLyBOVU1fRVhQUyksDQogICAgbGVuZ3RoLm91dCA9IG5yb3coZGF0YV9wZXJfZ3JvdXApDQogICkpDQoNCiMgc2V0IHVwIHBsb3QNCnAgPC0gZGF0YV9wZXJfZ3JvdXAgJT4lDQogIGdncGxvdCgNCiAgICBhZXMoDQogICAgICB4ID0gZHVtbXlfeCwgeSA9IG1lYW5fZGV2aWF0aW9uLCBjb2xvdXIgPSBleHBlcmltZW50DQogICAgKQ0KICApICsNCiAgdGhlbWVfY2xhc3NpYygpICsNCiAgIyB0aGVtZShsZWdlbmQucG9zaXRpb24gPSAibm9uZSIpICsNCiAgbGFicygNCiAgICB4ID0gIlRyaWFsIE51bWJlciIsDQogICAgeSA9ICJBYnNvbHV0ZSBUYXJnZXQgRXJyb3IgKGNtKSINCiAgKQ0KDQojIGFkZCBob3Jpem9udGFsIGxpbmVzDQpwIDwtIHAgKw0KICBnZW9tX2hsaW5lKA0KICAgIHlpbnRlcmNlcHQgPSBjKDAsIDQwKSwgbGluZXdpZHRoID0gMC40LA0KICAgIGNvbG91ciA9ICIjQ0NDQ0NDIiwgbGluZXR5cGUgPSAic29saWQiDQogICkgKw0KICBnZW9tX2hsaW5lKA0KICAgIHlpbnRlcmNlcHQgPSBjKDIwKSwgbGluZXdpZHRoID0gMC40LA0KICAgIGNvbG91ciA9ICIjQ0NDQ0NDIiwgbGluZXR5cGUgPSAiZGFzaGVkIg0KICApDQoNCiMgYWRkIGNvbmZpZGVuY2UgaW50ZXJ2YWxzIGFuZCBkYXRhIHBvaW50cw0KZm9yICh1bmlxdWVfdGVzdF90eXBlIGluIHVuaXF1ZShkYXRhX3Blcl9ncm91cCR0ZXN0X3R5cGUpKSB7DQogICMgZ2V0IHRoZSBkYXRhIGZvciB0aGlzIGJsb2NrDQogIHRvX3Bsb3RfZGF0YSA8LSBmaWx0ZXIoZGF0YV9wZXJfZ3JvdXAsIHRlc3RfdHlwZSA9PSB1bmlxdWVfdGVzdF90eXBlKQ0KDQogIHAgPC0gcCArIGdlb21fcmliYm9uKA0KICAgIGRhdGEgPSB0b19wbG90X2RhdGEsDQogICAgYWVzKA0KICAgICAgeW1pbiA9IG1lYW5fZGV2aWF0aW9uIC0gY2lfZGV2aWF0aW9uLA0KICAgICAgeW1heCA9IG1lYW5fZGV2aWF0aW9uICsgY2lfZGV2aWF0aW9uLA0KICAgICAgZmlsbCA9IGV4cGVyaW1lbnQNCiAgICApLCBjb2xvdXIgPSBOQSwgYWxwaGEgPSAwLjMNCiAgKSArIGdlb21fbGluZSgNCiAgICBkYXRhID0gdG9fcGxvdF9kYXRhDQogICkNCn0NCg0KIyBzZXQgY29sb3VyIHBhbGV0dGUNCnAgPC0gcCArDQogIHNjYWxlX2NvbG91cl9tYW51YWwodmFsdWVzID0gcGFsbGV0ZV9saXN0KSArDQogIHNjYWxlX2ZpbGxfbWFudWFsKHZhbHVlcyA9IHBhbGxldGVfbGlzdCkNCg0KZ2dwbG90bHkocCkNCmBgYA0KDQoNCg0KIyBBbmltYXRlZCBTdXJmYWNlIEZvbGxvdy11cA0KDQojIyMgUGxvdCBBTkdVTEFSIERFVklBVElPTlMgKGhhbmQgYW5nbGVzKQ0KYGBge3IsIGZpZy53aWR0aD0xMCwgZmlnLmhlaWdodD02LCBvdXQud2lkdGg9IjEwMCUifQ0KIyBpc29sYXRlIGFuaW1hdGVfc3VyZmFjZSBleHANCmRhdGFfcGVyX2dyb3VwIDwtIG9tbmlidXNfZGYgJT4lDQogIGZpbHRlcihleHBfbGFiZWwgPT0gImFuaW1hdGVfc3VyZmFjZSIpICU+JQ0KICBncm91cF9ieShwcmlvcl9hbmltLCBibG9ja19udW0sIHRyaWFsX251bV9pbl9ibG9jaywgdHJpYWxfbnVtKSAlPiUNCiAgc3VtbWFyaXNlKA0KICAgIG1lYW5fZGV2aWF0aW9uID0gbWVhbih0aHJvd19kZXZpYXRpb24pLA0KICAgIGNpX2RldmlhdGlvbiA9IHZlY3Rvcl9jb25maW50KHRocm93X2RldmlhdGlvbikNCiAgKQ0KDQojIG9yZGVyIHRoZSBmYWN0b3JzIGZvciBhc3NpZ25pbmcgY29sb3VyIHBhbGxldHMNCmRhdGFfcGVyX2dyb3VwJHByaW9yX2FuaW0gPC0gZmFjdG9yKA0KICBkYXRhX3Blcl9ncm91cCRwcmlvcl9hbmltLA0KICBsZXZlbHMgPSBjKA0KICAgICJub25lIiwgImhhbGZfYW5pbSIsICJmdWxsX2FuaW0iDQogICkNCikNCg0KIyBzZXQgdXAgcGxvdA0KcCA8LSBkYXRhX3Blcl9ncm91cCAlPiUNCiAgZ2dwbG90KA0KICAgIGFlcygNCiAgICAgIHggPSB0cmlhbF9udW0sIHkgPSBtZWFuX2RldmlhdGlvbiwNCiAgICAgIHltaW4gPSBtZWFuX2RldmlhdGlvbiAtIGNpX2RldmlhdGlvbiwNCiAgICAgIHltYXggPSBtZWFuX2RldmlhdGlvbiArIGNpX2RldmlhdGlvbg0KICAgICkNCiAgKSArDQogIHRoZW1lX2NsYXNzaWMoKSArDQogICMgdGhlbWUobGVnZW5kLnBvc2l0aW9uID0gIm5vbmUiKSArDQogIGxhYnMoDQogICAgeCA9ICJUcmlhbCBOdW1iZXIiLA0KICAgIHkgPSAiVGhyb3cgQW5nbGUgKMKwKSINCiAgKQ0KDQojIGFkZCBob3Jpem9udGFsIGxpbmVzDQpwIDwtIHAgKw0KICBnZW9tX2hsaW5lKA0KICAgIHlpbnRlcmNlcHQgPSBjKDAsIC0zMCksIGxpbmV3aWR0aCA9IDAuNCwNCiAgICBjb2xvdXIgPSAiI0NDQ0NDQyIsIGxpbmV0eXBlID0gInNvbGlkIg0KICApICsNCiAgZ2VvbV9obGluZSgNCiAgICB5aW50ZXJjZXB0ID0gYygtMTUpLCBsaW5ld2lkdGggPSAwLjQsDQogICAgY29sb3VyID0gIiNDQ0NDQ0MiLCBsaW5ldHlwZSA9ICJkYXNoZWQiDQogICkNCg0KIyBwIDwtIHAgKw0KIyAgIHNjYWxlX3lfY29udGludW91cygNCiMgICAgIGxpbWl0cyA9IGMoLTEwLCAzNSksDQojICAgICBicmVha3MgPSBjKDAsIDE1LCAzMCksDQojICAgICBsYWJlbHMgPSBjKDAsIDE1LCAzMCkNCiMgICApICsNCiMgICBzY2FsZV94X2NvbnRpbnVvdXMoDQojICAgICBsaW1pdHMgPSBjKDAsIDE4MCksDQojICAgICBicmVha3MgPSBjKDAsIDYwLCAxMjAsIDE4MCksDQojICAgICBsYWJlbHMgPSBjKDAsIDYwLCAxMjAsIDE4MCkNCiMgICApDQoNCiMgc2V0IGZvbnQgc2l6ZSB0byAxMQ0KcCA8LSBwICsNCiAgdGhlbWUodGV4dCA9IGVsZW1lbnRfdGV4dChzaXplID0gMTEpKQ0KDQojIHJlcGVhdCBmb3IgcHJpb3JfYW5pbSA9PSAiaGFsZiIsICJmdWxsIiBhbmQgIndhaXQiDQpmb3IgKHVuaXF1ZV9wcmlvcl9hbmltIGluIHVuaXF1ZShkYXRhX3Blcl9ncm91cCRwcmlvcl9hbmltKSkgew0KICAjIGdldCB0aGUgZGF0YSBmb3IgdGhpcyBibG9jaw0KICB0b19wbG90X2RhdGEgPC0gZmlsdGVyKGRhdGFfcGVyX2dyb3VwLCBwcmlvcl9hbmltID09IHVuaXF1ZV9wcmlvcl9hbmltKQ0KICAjIGxvb3AgdGhyb3VnaCB0aGUgdW5pcXVlIGJsb2NrcyBpbiB0b19wbG90X2RhdGENCiAgZm9yIChibG9jayBpbiB1bmlxdWUodG9fcGxvdF9kYXRhJGJsb2NrX251bSkpIHsNCiAgICAjIGdldCB0aGUgZGF0YSBmb3IgdGhpcyBibG9jaw0KICAgIGJsb2NrX2RhdGEgPC0gZmlsdGVyKHRvX3Bsb3RfZGF0YSwgYmxvY2tfbnVtID09IGJsb2NrKQ0KICAgICMgYWRkIHRoZSBkYXRhLCB1c2UgdGhlIHBhbGxldGVfbGlzdCB0byBnZXQgdGhlIGNvbG91cg0KICAgIHAgPC0gcCArIGdlb21fcmliYm9uKA0KICAgICAgZGF0YSA9IGJsb2NrX2RhdGEsDQogICAgICBhZXMoZmlsbCA9IHByaW9yX2FuaW0pLA0KICAgICAgY29sb3VyID0gTkEsIGFscGhhID0gMC4zDQogICAgKSArIGdlb21fbGluZSgNCiAgICAgIGRhdGEgPSBibG9ja19kYXRhLA0KICAgICAgYWVzKGNvbG91ciA9IHByaW9yX2FuaW0pDQogICAgKQ0KICB9DQp9DQoNCiMgc2V0IGNvbG91ciBwYWxldHRlDQpwIDwtIHAgKw0KICBzY2FsZV9jb2xvdXJfbWFudWFsKHZhbHVlcyA9IHBhbGxldGVfbGlzdCkgKw0KICBzY2FsZV9maWxsX21hbnVhbCh2YWx1ZXMgPSBwYWxsZXRlX2xpc3QpDQoNCiMgIyBzYXZlDQojIGlmIChzYXZlX3Bsb3RzKSB7DQojICAgZ2dzYXZlKA0KIyAgIHAsDQojICAgZmlsZW5hbWUgPSAiLi4vcGxvdHMvcGFwZXJfZmlncy9zcl8zMF90cmFpbmluZy5wZGYiLCBkZXZpY2UgPSAicGRmIiwNCiMgICBoZWlnaHQgPSA0LCB3aWR0aCA9IDYNCiMgICApDQojICAgfQ0KDQpnZ3Bsb3RseShwKQ0KIyBwDQpgYGANCg0KIyMjIFdhc2hvdXQgdHJpYWxzIGFmdGVyIGhhbGYgVlMgZnVsbCBhbmltYXRpb25zDQpgYGB7ciwgZmlnLndpZHRoPTEwLCBmaWcuaGVpZ2h0PTYsIG91dC53aWR0aD0iMTAwJSJ9DQojIGZpcnN0LCBpc29sYXRlIHRoZSBkYXRhDQpkYXRhIDwtIG9tbmlidXNfZGYgJT4lDQogIGZpbHRlcigNCiAgICBleHBfbGFiZWwgPT0gImFuaW1hdGVfc3VyZmFjZSIsDQogICAgYmFzZWxpbmVfYmxvY2sgPT0gRkFMU0UsDQogICAgdGVzdF90eXBlID09ICJ3YXNob3V0X2FuaW0iDQogICkNCg0KZGF0YV9wZXJfcHB0IDwtIGRhdGEgJT4lDQogIGdyb3VwX2J5KHBwaWQsIHByaW9yX2FuaW0sIHRyaWFsX251bV9pbl9ibG9jaykgJT4lDQogIHN1bW1hcmlzZSgNCiAgICBwcHRfbWVhbl9kZXZpYXRpb24gPSBtZWRpYW4odGhyb3dfZGV2aWF0aW9uKSwNCiAgICBwcHRfY2lfZGV2aWF0aW9uID0gdmVjdG9yX2NvbmZpbnQodGhyb3dfZGV2aWF0aW9uKSwNCiAgICBuID0gbigpDQogICkNCg0KZGF0YV9wZXJfZ3JvdXAgPC0gZGF0YV9wZXJfcHB0ICU+JQ0KICBncm91cF9ieShwcmlvcl9hbmltLCB0cmlhbF9udW1faW5fYmxvY2spICU+JQ0KICBzdW1tYXJpc2UoDQogICAgbWVhbl9kZXZpYXRpb24gPSBtZWFuKHBwdF9tZWFuX2RldmlhdGlvbiksDQogICAgY2lfZGV2aWF0aW9uID0gdmVjdG9yX2NvbmZpbnQocHB0X21lYW5fZGV2aWF0aW9uKSwNCiAgICBuID0gc3VtKG4pDQogICkNCg0KIyBzZXQgdXAgcGxvdA0KcCA8LSBkYXRhX3Blcl9ncm91cCAlPiUNCiAgZ2dwbG90KA0KICAgIGFlcygNCiAgICAgIHggPSB0cmlhbF9udW1faW5fYmxvY2ssIHkgPSBtZWFuX2RldmlhdGlvbiwNCiAgICAgIGNvbG91ciA9IHByaW9yX2FuaW0sIGZpbGwgPSBwcmlvcl9hbmltDQogICAgKQ0KICApICsNCiAgdGhlbWVfY2xhc3NpYygpICsNCiAgIyB0aGVtZShsZWdlbmQucG9zaXRpb24gPSAibm9uZSIpICsNCiAgbGFicygNCiAgICB4ID0gIlRyaWFsIE51bWJlciBpbiBCbG9jayIsDQogICAgeSA9ICJUaHJvdyBBbmdsZSAowrApIg0KICApDQoNCiMgYWRkIGhvcml6b250YWwgbGluZXMNCnAgPC0gcCArDQogIGdlb21faGxpbmUoDQogICAgeWludGVyY2VwdCA9IGMoMCwgLTMwKSwgbGluZXdpZHRoID0gMC40LA0KICAgIGNvbG91ciA9ICIjQ0NDQ0NDIiwgbGluZXR5cGUgPSAic29saWQiDQogICkgKw0KICBnZW9tX2hsaW5lKA0KICAgIHlpbnRlcmNlcHQgPSBjKC0xNSksIGxpbmV3aWR0aCA9IDAuNCwNCiAgICBjb2xvdXIgPSAiI0NDQ0NDQyIsIGxpbmV0eXBlID0gImRhc2hlZCINCiAgKQ0KDQojIGFkZCBkYXRhIHBvaW50cw0KcCA8LSBwICsgZ2VvbV9iZWVzd2FybSgNCiAgZGF0YSA9IGRhdGFfcGVyX3BwdCwNCiAgYWVzKA0KICAgIHkgPSBwcHRfbWVhbl9kZXZpYXRpb24sDQogICAgY29sb3VyID0gcHJpb3JfYW5pbQ0KICApLA0KICBzaXplID0gMSwgZG9kZ2Uud2lkdGggPSAwLjUNCikgKyBnZW9tX3JpYmJvbigNCiAgYWVzKA0KICAgIHltaW4gPSBtZWFuX2RldmlhdGlvbiAtIGNpX2RldmlhdGlvbiwNCiAgICB5bWF4ID0gbWVhbl9kZXZpYXRpb24gKyBjaV9kZXZpYXRpb24NCiAgKSwNCiAgY29sb3VyID0gTkEsIGFscGhhID0gMC4zDQopICsgZ2VvbV9saW5lKCkNCg0KDQojIHNldCBjb2xvdXIgcGFsZXR0ZQ0KcCA8LSBwICsNCiAgc2NhbGVfY29sb3VyX21hbnVhbCh2YWx1ZXMgPSBwYWxsZXRlX2xpc3QpICsNCiAgc2NhbGVfZmlsbF9tYW51YWwodmFsdWVzID0gcGFsbGV0ZV9saXN0KQ0KDQpnZ3Bsb3RseShwKQ0KYGBgDQoNCg0KYGBge3IsIGluY2x1ZGU9RkFMU0V9DQojIyBFTVBUWQ0KYGBgDQoNCg0KIyBEZXByZWNhdGVkIChEb24ndCBSdW4pDQoNCiMjIFN1Y2Nlc3MgbWFuaWZvbGRzDQojIyMgV2l0aG91dCBhbnkgdGlsdHMNCmBgYHtyLCBmaWcud2lkdGg9OSwgZmlnLmhlaWdodD0yMH0NCiMgZ2dwbG90bHkocGxvdF9zdWNjZXNzX21hbmlmb2xkX25vX3RpbHQoKSkNCnBsb3Rfc3VjY2Vzc19tYW5pZm9sZF9ub190aWx0KCkNCmBgYA0KDQojIyMgV2l0aCB0aWx0IHByZXNlbnQNCmBgYHtyLCBmaWcud2lkdGg9OSwgZmlnLmhlaWdodD0yMH0NCmdncGxvdGx5KHBsb3Rfc3VjY2Vzc19tYW5pZm9sZF90aWx0KCkpDQpgYGANCg0KDQpgYGB7ciwgaW5jbHVkZT1GQUxTRX0NCk5VTEwNCmBgYA0KDQo=